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electrode  separation. 

In  order  to  understand  hole  transport  and  mobility  in  SiO  ,  we.  have 
studied  optical  transitions  by  using  uv  reflectivity  to  investigate 
structures  near  the  band  edge  as  a  function  of  the  Si-O-Si  bond  angle. 
Structure  is  seen  in  the  optical  data  of  cristoballite  which  indicates 
that  the  uppermost  narrow  valence  band  has  a  finite  bandwidth  due  to 
wavef unction  overlap. 

Phosphorus  precipitates  were  found  on  the  interface  between  P-doped 
silicon  and  thermally  grown  SiO„.  The  technique  of  scanning  internal 
photoemission  was  used  to  form  an  image  of  the  phosphorus  islands  whose 
size  and  spacing  depend  upon  the  areal  density  of  phosphorus  in  the 
oxidized  silicon.  The  islands  were  found  to  completely  cover  the 
silicon  surface  for  an  areal  density  of  5  x  lO^/cm. 

Resistive  switching  was  studied  in  Ge1  ,.Tefl[.  by  transmission  electron- 
microscopy.  Three  resistive  states  were  observed,  consisting  of  a  wholly 
amorphous  structure,  Te  crystallites  in  an  amorphous  GeTe  matrix,  and  Te 
crystallites  in  a  crystalline  GeTe  matrix.  Switching  was  observed 
between  the  latter  two  but  structural  degradation  resulted. 

In  the  area  involving  the  structural  instability  of  glass-metal 
interfaces,  we  have  performed  a  systematic  study  of  interfacial  reactions 
jetween  thermally-grown  S1C>2  and  three  transition  metals;  V,  Tis  and  Nb. 

We  found  that  there  interfaces  became  unstable  at  temperatures  above 
700  C  and  destruction  of  the  structure  of  SiO  followed.  This  is 
significant  because  active  transition  metals  are  often  used  to  increase 
the  adhesion  of  a  second  metal  to  the  surface  of  SiO?. 
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THE  PHYSICS  OF  INTERFACE  INTERACTIONS  RELATED  TO 
RELIABILITY  OF  FUTURE  ELECTRONIC  DEVICES 

INTRODUCTION 

The  overall  purpose  of  our  work  continues  to  be  the  investigation 
of  fundamental  phenomena  that  directly  impact  the  performance  of 
proiected  electronic  device  configurations.  We  have  been  probing  and 
reporting  on  various  possible  approaches  to  understanding  the  many 
phenomena  which  can  influence  the  reliability  of  integrated  devices 
including  such  areas  as  polarization  effects  on  interface  contact 
barriers  and  carrier  injection  into  Si02;  dielectric  breakdown  mechan¬ 
isms  and  models,  band  structure  in  insulation;  radiation-indued 
charge  in  insulators;  trapping  states  in  insulation  interfaces  such 
as  MHOS  structures;  switching  and  conduction  mechanisms  in  doped 
transition  metal  oxides;  drift  of  charged  carriers  through  Si02  and 
impurity  Induced  crystallization;  and  reactions  at  metal-Si  and  metal- 
Si02  interfaces. 

For  the  near  future,  we  will  concentrate  on  those  programs  which 
have  proven  to  be  of  fundamental  scientific  importance  and  which  under¬ 
lie  the  significant  degradation  modes  that  can  be  foreseen.  This 
provides  the  best  chance  for  future  extension  of  new  device  concepts. 
These  include  a  more  concerted  effort  theoretically  to  describe  inter¬ 
face  polarization  effects  and  contact  barrier  drift,  continued  expan¬ 
sion  of  the  theoretical  dielectric  breakdown  model  and  its  implications, 
description  of  radiation-induced  charge  carriers  in  Si02,  and  reaction 
kinetics  at  metal-Si02  (or  Si)  interfaces. 
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O  We  have  developed  an  understanding  of  precisely  what  happens  in 
the  initial  stages  of  dielectric  breakdown  in  Si(>2.  A  simple  model 
for  dielectric  breakdown  was  formulated  on  the  basis  of  electron 
inlection  from  blocking  contacts  and  subsequent  impact  ionization  of 
the  lattice.  Several  predictions  can  be  made,  including  the  increase 
average  breakdown  field  found  for  ultra-thin  films  of  Si02>  Ionizing 
radiation  is  predicted  to  reduce  drastically  the  breakdown  strength 

of  Si02. 

O  The  breakdown  model  predicts  ^n  increase  in  the  breakdown  field 
of  greater  than  100%  for  SiO  films  100  X  thick.  For  such  thin  films, 
the  negative  reistance  portion  of  the  curve  disappears.  This  means 
that  thin  films  about  100  R  thick,  would  have  superior  qualities  as  a 
gate  insulator. 

O  For  n-type  silicon  or  for  n-channel  MOS  devices,  the  p.iosphorus 
dopant  from  the  silicon  was  found  to  precipitate  onto  any  silicon 
surface  which  has  been  oxidized.  This  surface  layer  or  phosphorus 
can  cause  a  leakage  channel  to  appear  near  the  surfaces  of  heavily 
phosphorus  doped  silicon.  As  a  result,  the  surface  leakage  path 
can  destroy  the  operation  of  an  ultra-small  MOS  transistor  because 
the  leakage  path  bridges  the  source  and  drain  regions.  It  is 
expected  that  arsenic,  the  other  n-type  dopant,  displays  the  same 
behavior. 

O  We  have  found  that  the  optical  properties  of  SiO  do  not  depend 
significantly  upon  the  alio  tropic  form  of  the  material.  The  res^lt:s 
indicate  a  valence  band  width  of  about  1  eV  for  the  top  valence  an 
of  Si02  crystoballite. 

O  We  have  found  that  certain  metals  can  cause  reliability  problems 
when  used  as  metallization  on  SiO  glass,  especially  where  the 
combination  is  subjected  to  elevated  temperatures.  Metals  that  form 
both  stable  silicides  and  oxides,  such  as  V,  Ti,  Cr,  Nb,  and  Ta,  may 
decompose  SiO  thin  films  over  a  period  of  time,  and  thereby  destroy 
the  dielectric  properties  of  the  insulating  film. 

O  In  the  area  of  Schottky  bs-riers,  we  have  finished  a  review  of 
our  current  understanding  of  sij icide  formation  and  stability.  We 
found  that  there  is  a  distinct  difference  between  the  silicides  formed 
by  the  metal-silicon  reaction  and  by  metal-  glass  reaction;  the  latter 
produce  a  metal-rich  phase  while  the  former  produce  a  silicon-rich 
phase.  This  shows  the  signficance  of  having  a  thin  oxide  layer  on 
the  contact  reaction  between  a  Si  wafer  and  a  metallic  films. 
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I.  THE  RELIABILITY  OF  SEMICONDUCTOR- INSULATOR  INTERFACES 

We  have  completed  an  investigation  of  a  physical  model  of  dielectric 
breakdown  in  Si02  thin  films,  including  a  study  of  the  influence  of  the 
adjustable  parameters  and  a  study  of  several  important  predictions. 

This  is  the  first  quantitative  model  to  describe  exactly  what  occurs 
during  the  critical  initial  stages  of  dielectric  breakdown  in  the 
dielectrics  used  in  integrated  circuits.  The  negative  resistance  type 
of  breakdown  is  due  to  a  distortion  of  the  electric  field  caused  by 
positive  charge  left  in  the  insulator  by  the  small  amount  of  impact 
ionization  produced  by  hot  electrons.  From  this  study,  we  find  that 
the  only  parameter  which  has  a  significant  effect  on  the  calculated 
breakdown  strength  is  the  electron-phonon  scattering  length.  A  good 
fit  to  the  experimental  breakdown  data  for  films  below  1000  X  was 
obtained  for  a  scattering  length  of  1.34  X.  The  analysis  was  expanded 
to  include  relative  effects  of  the  mean  free  path  for  impact  ion¬ 
ization  of  the  valence  band  edge  as  well  as  the  recombination  cross 
section  a  as  calculated  by  several  models.  Neither  were  as  important 
In  determining  breakdown  as  the  electron-phonon  scattering  length. 

The  theory  predicts  the  experimentally  observed  increase  in  break¬ 
down  strength  for  films  thinner  than  ^  500  X.  The  other  important 
predictions  need  experimental  verification:  the  threshold  current  for 
breakdown  is  found  to  increase  dramatically  with  a  decreasing  film 
thickness,  and  the  intrinsic  breakdown  voltage  decreases  in  the  presence 
of  ionizing  radiation. 

By  scanning  internal  photoemission  measurements,  we  have  seen 
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evidence  for  the  precipitation  of  phosphorus  on  Si-Si02  interfaces 
grown  on  phosphorus  doped  silicon.  Apparently,  the  phosphorus  that 
dissolves  in  the  Si02  during  oxidation  of  the  doped  silicon  precipitates 
onto  the  interface  during  the  period  when  the  silicon  is  cooled  from 
the  oxidation  temperature.  Islands  of  phosphorus  were  seen  as  dark 
spots  (low  photoyield)  on  photoemission  images  of  sodium  stained  Si-Si02 
interfaces.  These  spots  are  absent  from  interfaces  grown  on  lightly 
phosphorus  doped  silicon  2  x  1018/cm3> ,  while  they  completely  cover 
an  interface  grown  on  a  very  heavily  doped  wafer  (n,  5  x  1020/cm3) . 
Samples  in  the  intermediate  range  of  doping  show  phosphorus  rich  islands 
on  the  Si-Si02  interface,  with  a  diameter  of  up  to  25  microns.  The 
significant  fact  is  that  silicon  surfaces  are  degraded  by  phosphorus 
which  is  introduced  during  the  thermal  oxidation  of  heavily  doped 
silicon. 


A.  IMPACT  IONIZATION  MODEL  FOR  DIELECTRIC 


INSTABILITY  AND  BREAKDOWN 


T.  H.  D i  Stefano  and  M.  Shatzkes 

Dielectric  breakdown  ir.  wide-bandgap  insulators  is  an  old  problem 
which  has  long  been  considered  in  quite  general  terms  from  several  points 
of  view.  These  treatments  determine  criteria  for  electronic  instabil¬ 
ity  and  dielectric  breakdown  in  terms  of  the  behavior  of  a  single  electron 
in  an  infinite  insulator  which  supports  a  uniform  electric  field.  More 
recent  theoretical  work  has  resulted  in  a  few  microscopic  models  for  diel¬ 
ectric  breakdown  in  finite  systems.  The  models  fall  into  two  classifica¬ 
tions:  an  instability  due  to  thermal  runaway3  in  a  system  containing  a 
relatively  low  energy  excited  electronic  state;  or  an  electronic  avalanche4’5 
type  of  dielectric  breakdown  in  which  hot  electrons  impact  bound  electrons 
and  ionize  them  from  the  valence  band.  By  either  a  thermal  or  an  avalanche 
model,  the  current-voltage  relationship  is  found  to  have  a  negative  resis¬ 
tance  or  dielectric  instability  due  to  a  non-local  regenerative  mechanism. 

The  avalanche  model  for  electronic  breakdown  developed  by  O'dwyer4  employs 
a  distortion  of  the  electric  field  in  the  insulator  to  obtain  the  non-local 
feedback  necessary  for  the  instability.  In  effect,  the  positive  charge 
carriers  produced  during  an  avalanche  drift  to  the  cathode  where  they  en¬ 
hance  the  electric  field  which,  in  turn,  leads  to  an  increased  electron 
injection  current  and  a  greater  rate  of  impact  ionization.  The  O'dwyer 
model  is  based  on  the  assumption  that  the  holes  produced  by  impact  ioniza¬ 
tion  remain  relatively  mobile  and  that  the  electron  energy  distribution 
depends  only  upon  the  local  electric  field  but  not  upon  distance  from  the 
cathode.  This  set  of  assumptions  is  valid  for  many  semiconductors,  but  it 
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does  not  apply  well  t0  thin  fil„3  of  those  ulde.bandgap  lnsulator3  ^  __ 
very  low  hole  mobility. 

We^ present  a  simple  theory  for  dlslectric  breakdown  in  materials  such 
as  S1026  for  which  the  hole  mobility  is  relatively  small.  The  calculations 
are  based  on  impact  Ionisation  of  the  insulator  lattice  by  electrons  injec¬ 
ted  from  the  cathode  and  subsequently  heated  by  the  electric  field.7  The 
theory  is  developed  lor  the  idealised  case  in  which  hole  motion  is  negli¬ 
gible  during  the  rapid  breakdown  process.  It  is  assumed  that  during  the 
short  time  required  for  breakdown,  hmes  are  cleared  from  the  Insulator 
largely  by  recombination  and  that  drift  in  the  electric  field  does  not  ser¬ 
iously  influence  the  hole  distribution  or  density  in  the  Insulator.  As  a 
specific  example  of  this  model,  the  breakdown  voltage  and  the  current- 
voltage  relationship  were  calculated  for  the  case  of  a  thin  layer  of  SiOj. 
The  breakdown  voltage  was  found  to  depend  significantly  upon  only  one  of2 
the  adjustable  parameters,  the  one  dimensional  energy  loss  length.  The 
ret ombination  cross  section  and  the  ionisafon  cross  section  have  a  rather 
small  influence  on  the  calculated  breakdown  voltage.  The  dependence  of  the 

average  electric  field  at  breakdown  on  thickness  is  calculated  and  found  to 
increase  rapidl  belov?  200A. 

The  principal  features  of  the  model  are  outlined  in  Fig.  1,  which 
shows  the  valence  and  conduction  bands  in  the  insulator  during  the  break¬ 
down  process.  Since  the  bandgap8  of  Si02  is  9.0  eV,  the  threshold  for 
impact  ionization  is  9  eV  above  the  bottom  of  the  conduction  band,  as 
marked  by  the  dashed  line  in  the  ,igure.  Electrons  are  injected  into  the 
insulator  by  Fowler-Nordheim  tunneling  from  one  of  the  electrodes.  The 
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tunneling  current  density  J  is  of  the  form9, 

J  -  (q3E2/16Tr2-fi2$)e"t(32m*)1/2  $3/2/3«eE] 

wh.,e  .  is  the  interface  energy  barrier,  m*  is  the  effective  mass  in  the 
insuiator,  and  E  la  the  electric  field  at  the  cathode.  For  SiO,, ,  „*-.41 
and  *-4.2  eV.  As  will  be  shown,  the  results  of  the  calculation  of  the 
field  distortion  in  the  insulator  are  independent  of  the  level  of  injected 
.current.  Since  the  solution  is  insensitive  to  the  level  of  injected  cur¬ 
rent  for  hi^h  »  contacts  ami  low  tesperatures ,  we  will  use  the  approximate 

form  for  the  Fowler-Nordheim  current,  and  neglect  temperature  and  image 
force  lowering  corrections. 

The  injected  electrons  are  accelerated  by  the  electric  field  so  that 
after  traversing  a  distance  x,  a  high  energy  tail  of  the  electron  energy 
distribution  has  sufficient  energy  to  ionize  the  lattice.  Electrons  above 
the  ionization  threshold  are  assumed  to  produce  impact  ionization  with  a 
mean  length  Each  impact  ionization  event  leaves  behind  a  low-mobility 

valence  band  hole.  In  this  model,  it  is  assumed  that  the  hole  moves  an 
insignificant  distance  before  it  is  anihilated  by  recombination  with  another 
injected  electron.  As  a  consequence  of  repeated  impact  ionization,  a  re¬ 
sidual  cloud  of  positive  space  charge  is  built  up  and  the  potential  is 
distorted  so  as  to  enhance  the  cathode  field.  In  turn,  the  enhanced  field 
leads  to  impact  ionization  nearer  the  cathode,  moving  the  sheet  of  positive 
charge  closer  to  the  cathode.  The  holes  downstream  from  the  major  ioniza¬ 
tion  region  are  removed  by  recombination.  By  a  regeneration  process,  the 
positive  charge  sheet  moves  closer  to  the  cathode,  leading  to  an  increase 


7 


in  the  cathode  field  and  a  resulting  increase  in  the  injected  electron 
current.  At  a  critical  point,  an  increase  in  current  is  accompanied  by  a 
decrease  in  total  potential  drop  across  the  thin  film.  This  point,  at  the 
beginning  of  the  negative  resistance  regime,  marks  the  onset  of  a  dielec¬ 
tric  breakdown  instability. 

In  order  to  calculate  the  current-voltage  relationship  and  to  find 
the  dielectric  breakdown  voltage  of  a  thin  dielectric,  the  rate  of  impact 
ionization  is  first  determined.  The  problem  is  complicated  by  the  fact 
that  the  electron  energy  distribution  does  not  reach  an  equilibrium  form 
in  the  sharply  varying  field  of  the  insulator.  It  is  both  possible  and 
necessary  to  determine  the  spacial  evolution  of  the  electron  distribution, 
with  the  aid  of  a  few  simplifying  approximations.  The  electron  scattering 
events  are  assumed  to  be  independent  and  to  involve  the  loss  of  one  LO 

phonon  per  collision  which  is  assumed  to  be  predominantly  e  forward  scat¬ 

tering  event.  The  electron  energy  distribution  is  obtained  by  considering 
the  one  dimensional  projection  of  the  scattering  events  to  be  a  Poisson 
process.  That  is,  the  electron  moves  an  average  distance  A  in  the  x 
direction  between  independent  scattering  events.  The  parameter  A  is  the 

one  dimensional  projection  of  the  phonon  emission  length.  This  A  is  not 

to  be  confused  with  the  mean  free  path  or  momentum  randomization  length; 
indeed,  the  two  can  be  quite  different  for  the  case  of  a  predominantly  for¬ 
ward  scattering  event.  The  probability  that  an  electron  will  suffer  n 
collisions  in  traversing  a  distance  A  is, 

P(n,x)  =  I(£)Vx/X.  (2) 
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A  discrete  energy  distribution  curve  at  point  x  is  obtained  by  allowing 
the  electron  to  lose  one  phonon  or  Ku)  in  each  collision.  At  x,  the  prob¬ 
ability  P  is  a  function  of  the  energy  E, 

p(x,E)  =  (-e/KuYT  e  ’  (3) 

or,  by  the  central  limit  theorem,  the  energy  distribution  is 


D(x,E) 


~  (l/Km) 

/2tt(x/X) 


(E/Kui  +  x/A)^ 
e  2(x/X) 


(4) 


where  E  is  the  energy  with  respect  to  the  initial  energy. 

For  the  case  in  which  an  electron  ionizes  the  lattice  as  soon  as  it 

attains  an  energy  of  above  the  conduction  band  edge,  the  rate  of  impact 

ionization  or  production  of  positive  charge  is  then, 

o 


dx 


/ 

$(x)  +  E 


D(x,E)  dE. 


(5) 
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In  the  equilibrium  condition,  the  rate  of  production  of  positive  charge  by 
ionization  is  equal  to  that  of  destruction  by  recombination. 


+  ■  e 
r 


^2+  ^  a (E)  D(X,E)  dE 


(6) 


where  o(E)  is  the  energy  dependent  recombination  cross  section.  The  hole 
charge  density  p_^  is^found  to  b'i  independent  of  injected  current  J, 


d  ^  D(.XjE)  dE 

p  ■  dx  <J>(x)  +  E 

_E _ _ 


,  (7) 


,  ♦ 00  o 

o(E  )  /  D(:.,E)  dE  +  ;  a{E)  D(x,E)  dE 

-oo  <f>00 

where  o(E)  is  the  energy  dependent  recombir at  Ion  ~ross  section  and  ET  is 
the  limiting  energy  of  electrons  near  the  conduction  band  bottom.  Here, 
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the  tail  of  the  electron  distribution  which  is  calculated  to  be  at  or  be¬ 
low  the  conduction  band  bottom  is  assumed  to  have  a  limiting  energy10  of 

V4R“/3-  Now  the  potential  p(x)  can  be  found  from  a  double  integration 

of  the  net  charge  in  the  insulator. 

<p(x)  =  j  /  dx  /  dx  (P+  -  p_).  (8) 

In  the  usual  situation,  p_  is  negligible  in  comparison  to  p+.  The  po¬ 
tential  V  drop  across  a  film  of  thickness  d  is  found  by  the  integration 

the  current  corresponding  to  the  V(d)  is  obtained  from  the  Fowler- 
Nordheim  equation(1)  evaluated  for  the  electric  field  at  the  cathode.  From 

these  two  values  for  Ax=d,  the  J-V  relationship  is  calculated  for  an  insu¬ 
lator  of  thickness  d. 

Our  impact  ionization  model  is  applied  to  the  npecific  case  of  Si02> 
which  may  be  somewhat  less  than  ideal  for  two  reasons.  First,  although  the 
hole  mobility  in  Si02  is  quite  small  and  unknown  it  must  be  finite.  Since 
this  mobility  is  undetermined  at  present,  we  leave  it  out  of  the  calcula¬ 
tion;  however,  more  exact  solution  including  the  effect  of  hold  motion  are 
not  expected  to  influence  the  general  characteristics  of  the  model. 

Another  assumption  which  is  not  exact  for  Si02  is  the  phonon  energy  loss, 
which  is  taken  to  be  .the  LO  phonon  energy  Kiu=.153  eV.  At  electron  ener¬ 
gies  near  the  ionization  threshold,  phonons  from  other  parts  of  the 
Br illouin  zone  contribute  to  the  scattering  but  introduce  an  energy  loss 
less  than  the  full  .153  eV.  Since  the  relative  cross  sections  are  not 
known  for  scattering  by  the  various  types  of  phonons,  we  must  retain  the 
simple  Einstein  model  for  the  phonon  spectrum. 
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Milne-model  by  about  a  factor  of  106.  Yet  the  same  breakdown  voltage  is 

O  Q 

obtained  if  a  A=1.68  A  is  used  in  place  of  the  X=1.34  A  required  in  the 
Milne  model. 

In  figure  4,  the  average  field  at  breakdown,  the  breakdown  voltage 
divided  by  the  oxide  film  thickness,  is  plotted  as  a  function  of  oxide 
thickness  for  the  Milne  recombination  cross-section  and  1=1.34  A.  There 
is  little  change  in  the  average  field  for  the  thicker  films,  while  for 
films  thinner  than  ~  200  A  there  is  a  sharp  increase  in  the  average  field. 

In  this  calculations  we  have  retained  what  we  considered  to  be  the 
essential  characteristics  of  insulators  such  as  SiO,,  i.e.  the  relative 
immobility  of  the  hole,  and  a  single  active  optical  phonon.  Neglect  cf 
hole  motion  and  other  energy  loss  mechanisms  are  approximations  that  re¬ 
duced  the  complexity  of  the  problem,  but  may  have  introduced  some 
inaccuracy. 
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Figure  1:  Electrons  injected  from  the  cathode  lose  energy  via  phonon 
collisions,  so  that  the  electron  energy  distribution  at  a  position  within 
the  insulator  will  be  as  shown.  Those  electrons  with  kinetic  energy  larger 
than  the  ionization  energy  unde  go  impact  ionization,  losing  energy  in  the 
process;  an  electron  is  excited  into  the  conduction  band  from  the  valence 
band  leaving  an  immobile  hole.  Electrons  recombine  with  the  holes.  The 
positive  space  charge  causes  an  enhancement  of  the  field  at  the  cathode. 


14 


Figure  2:  Current  density  vs  applied  voltage  for  a  1000  %  Si02  layer.  For 
current  densitieu  below  the  onset  of  negative  resistance,  the  curve  is 
indeoendent  of  the  value  chosen  for  the  ionization  mean  free  path  X^.  In 
the  negative  resistance  regime  variations  with  X^  are  significant.  The 
oreakdown  voltage  is  the  maximum  applied  voltage  such  that  a  steady  state 
solution  for  the  current  is  possible. 
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Figure  3:  Current  density  vs  applied  voltage  for  two  models  °1 
hole  recombination.  In  both  cases  1  -  0  for  the  Milne  model  -=hjtf«™n 

mean  free  path  X  -  1.34  X,  while  for  the  Rose  model  X  -  1.68  A.  With  these 
choices  for  X,  the  breakdown  voltage  is  the  same  for  both  mgdels,  though 
recombination  cross  sections  differ  by  about  a  factor  of  .  , 

cross  section  being  smaller.  This  illustrates  the  extreme  insensitivity 
of  the  calculation  to  variations  in  the  recombination  cross  section. 
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B.  PHOSPHORUS  PRECIPITATION  ON  Si/Si02  INTERFACES 
T.  H.  DiStefano 

Scanning  internal  photoemission  has  been  used  to  study  the 
recently  discovered  phenomenon1 *2  that  phosphorus  is  found  in  a  thin 
layer  on  the  oxide  side  of  a  Si-SiO?  interface  formed  by  the  oxidation 
of  phosphous  doped  silicon.  It  is  well  knewn  that  the  silicon  surface 
of  each  interface  has  an  ennaneed  phosphorus  concentration.  The 
phosphorus  found  on  the  interface  is  different  in  that  it  comprises 
a  large  fraction  of  a  monolayer  of  phosphorus  in  a  layer  which  is 
10  8  or  less  in  thickness.  By  using  scanning  internal  photoemission, 
we  have  been  able  to  examine  the  distribution  of  phosphorus  on  the 
interface.  We  have  found  islands  of  phosphorus,  from  1  to  10  microns 
in  diameter  or  even  larger  in  some  cases,  distribution  across  the 

interface.  The  island  size  depends  upon  the  concentration  of  phosphorus 
in  the  substrate. 

The  interfaces  which  were  studied  were  prepared  by  thermal 
oxidation  of  the  (100)  surface  of  phosphorus  doped  silicon.  The 
oxidation  was  performed  by  the  standard  drv-wet-dry  technique  at 
1100°C,  to  obtain  a  film  of  Si02  to  a  thickness  of  1000  8.  Before 
oxidation,  the  silicon  surface  was  phosphorus  doped  by  diffusion  to  a 
depth  of  %  2  microns;  the  surface  of  the  silicon  was  cleaned  after 
the  phosphorus  diffusion  by  a  surface  etch  procedure.  The  photo¬ 
emission  measurements  were  performed  on  a  series  of  interfaces  grown 
on  silicon  with  a  range  of  phosphous  doping  levels. 

The  photoemission  images  are  shown  in  Figs.  1-4,  where  the  field 
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of  view  is  250  microns  by  250  microns.  Interfaces  grown  on  a  silicon 

19  2 

with  a  doping  level  of  about  5  x  10  P/cm  ,  shown  in  Figs.  1-2,  show 

phosphous  islands  about  10  microns  in  diameter.  For  a  substrate  doping 

level  of  about  2  x  10  P.cm  ,  the  phosphous  islands  are  much  smaller 

and  further  apart  as  seen  in  Figs.  3-4.  On  interfaces  grown  on  less 

heavily  doped  silicon  show  few  of  these  phosphous  islands,  while  those 
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grown  on  silicon  doped  to  lo  P/cm  are  completely  covered  with  the 
phosphorus  precipitate.  The  influence  of  the  phosphorus  doping  level 
is  summarized  in  Fig.  5,  which  shows  the  photocurrent  produced  by 
light  at  3250  X  for  interfaces  grown  on  silicon  with  a  range  of 
phosphorus  doping.  It  is  seen  that  the  average  photocurrent  drops 
rapidly  with  increasing  phosphorus  up  to  a  level  of  10  P/cm  ,  at 
which  point  the  photocurrer.t  disappears  completely.  This  phosphorus 
level  corresponds  to  a  concentration  of  5  x  10~^  P/cm2  on  the  inter¬ 
face,  if  it  is  assumed  that  all  of  the  phosphorus  in  the  oxidized 
layer  of  silicon  has  segregated  at  the  Si-Si02  interface. 
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Figure  1  and  2:  Scanning  Internal  Photoemission  maps  of  an  Si-SiC>2 
interface  gronw.on  silicon  (100)  with  a  Surface  phosphorus  doping 
level  of  5  x  10  P/cm  . 
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Figure  3  and  4:  Scanning  Internal  Photoemission  maps  of  a  Si-SiC^ 
interface  grown.,  on  silicon  (100)  with  a  surface  phosphorus  doping 
level  of  2  x  10  P/cro  . 
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Figure  5:  The  relative  photocurrent  measured  from  Si-SiO„  interfaces 
grown  on  (100)  silicon  with  a  range  of  phosphorus  doping  levels. 


II.  BAND  STRUCTURE  AND  SWITCHING  IN  INSULATORS 

From  our  theoretical  work  on  dielectric  breakdown  mechanisms, 
it  has  become  apparent  that  more  information  is  necessary  to  under- 
stand  the  role  of  holes  in  the  SiC^  formed  by  impact  ionization  or  by 
radiation.  Accumulation  of  positive  charge  causes  increased  field 
at  the  cathode  and  increased  carrier  injection.  The  stability  of  the 
holes  in  terms  of  their  mobility  is  key  to  the  predictions  of  the  theory. 
To  help  gain  needed  insight  into  hole  mobility,  it  is  desireable  to 
ascertain  the  valence  band  width.  To  do  this,  the  effects  of  crystal 
field  splitting  must  be  minimized.  This  was  accomplished  by  measuring 
the  reflectivity  of  the  allotropes  of  SiC^  characterized  by  different 
angles  of  the  Si-O-Si  bond:  a  quartz,  amorphous  quartz,  cristobalite. 

The  latter  case  should  not  show  crystal  field  splitting.  The  optical 
transitions  giving  rise  to  the  reflectivity  results  indicated  a 
valence  band  width  of  about  1  eV.  This  finite  width  implies  that  hole 
conduction  is  possible,  but  the  small  width  means  that  any  lattice 
defect,  such  as  unbridged  oxygen,  will  induce  trapping  of  any  holes 
in  the  valence  band  (Anderson  localization). 

The  second  paper  deals  with  change  in  resistance  state  of  Ge1cTeoc 
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films  by  switching  between  three  different  structures.  The  highest 
resistance  state  is  found  for  the  completely  amorphous  film.  When 
switching  to  a  lower  resistance  state  crystallites  of  Te  are  formed 
in  a  still  amorphous  matrix.  The  lowest  reistance  state  involves  a 
second  switch  in  which  the  matrix  crystallizes  to  GeTe.  Switching 
between  the  lowest  two  States  is  readily  accomplished.  However,  by 
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simultaneously  observing  the  structural  changes  occurring  during  the 
switching  process  by  transmission  electron  microscopy  it  can  be  con¬ 
cluded  that  permanent  degradation  of  the  film  increases  until  failure 
ultimately  takes  place.  This  seems  typical  of  other  switchable 
materials  we  have  investigated  in  the  past  such  as  NbjO,.  where  atomic 
rearrangement  accompanies  the  switching  process. 
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A.  OPTICAL  PROPERTIES  0 Tf  ALLOTROPIC  FORMS  OF  Si02 
T.  H.  DiStefano  and  P.  K.  Roy 

Until  recently,  little  has  been  reported  about  the  electronic 
structure  of  Si02  which  determines  several  physical,  optical  and 
dielectric  properties.  Knowledge  of  the  valence  band  (VB)  structure 
is  important  to  current  studies  of  radiation  damage,  electronic 
structure  of  various  defects,  and  properties  of  Si  (or  metal)-Si02 
interfaces.  Optical  properties  (ultraviolet  reflectance  spectra  from 
8  to  1.3  eV  in  the  present  case)  for  different  allotropes  of  Si02  will 
give  - 

(a)  detailed  information  about  the  structure  of  the  relatively 
unknown  VB  edge. 

(b)  variations  of  optical  and  dielectric  properties  as  a  function 
of  Si-O-Si  bond  angle,  0.  As  we  know  0  varies  from  120°  (coesite)  to 
180°  (for  cristobalite  or  tridymite)  with  crystalline  and  amorphous 
quartz  having  bond  angles  of  144°  and  150°,  respectively. 

(c)  The  width  of  the  top  of  the  VB  (AE,  Fig.  2),  which  plays 

a  very  important  role  in  determining  the  mobility  of  the  holes  created 
when  Si02  is  subjected  to  an  ionizing  radiation. 

The  structure  of  amorphous  Si02  era  be  understood  in  terms  of 
8-crystobalite,  in  which  the  SiO^  tetrahedra  occupying  the  points  of 
the  diamond  lattice  are  connected  by  straight  Si-O-Si  chains.  Oxygen 
in  this  lattice  occupies  the  sites  of  D3d  symmetry  shown  in  Fig.  1 
(a).  In  the  amorphous  Si02,  the  Si-O-Si  chain  is  bent  to  an  angle  of 
about  150°.  The  electronic  levels  in  the  amorphous  structure  are 
conveniently  described  with  the  0  2S  2p  orbitals  as  a  basis.  As 
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displayed  in  Fig.  1(b),  the  s  and  p  orbitals  of  the  0  ion  in  free 

space  are  separated  by  14.9  eV.  In  the  D3d  symmetry  of  B-crystobalite, 

the  2p6  degeneracy  is  split  to  produce  two  nonbonding  orbitals  (r&) 

which  are  perpendicular  to  the  Si-O-Si  chain,  and  one  bonding  orbital 

(r3)  which  links  the  Si  and  0  centers.  The  s  orbitals  remain  as  the 

T  level  in  D  .  In  going  to  the  amorphous  structure,  the  electronic 

levels  found  in  D„,  broaden  into  bands  as  shown  in  Fig.  1.  Also,  by 
Jd 

symmetry  considerations,  the  IV  level  is  shown  to  split  into  two  non¬ 
bonding  bands,  corresponding  to  nonbonding  orbitals  which  lie  perpen¬ 
dicular  to  and  parallel  to  a  plane  determined  by  the  Si-O-Si  chain. 

The  two  nonbonding  bands,  the  bonding  bands,  and  the  deep-lying  s  band 
described  in  terms  of  an  0"-ion  basis,  are  equivalent  to  those  found 
theoretically  by  Reilly  who  used  hybrid  Si  orbitals  and  atomic  0 
orbitals  as  a  basis.  Reilly's  determination  of  the  energy  levels  of 
SiC^  using  linear  combinations  of  atomic  and  molecular  orbitals  shows 
two  nonbonding  levels  at  the  valence-band  edge,  split  by  about  2  eV, 
with  a  bonding  level  at  about  12  eV  below  the  valence-band  edge.  The 
oxygen  2s  level  was  placed  at  about  24  eV  below  the  band  edge. 

Details  of  the  emission  spectra  from  SiC>2  can  be  identified  in 
terms  of  bands  in  Fig.  1(b).  The  two  peaks  at  the  valence-band  edge 
have  a  total  width  at  half-maximum  of  only  3.2  eV,  presumably  because 
of  the  small  overlap  between  the  nonbonding  orbitals.  The  bonding 
band,  seen  to  extend  from  about  -5.5  to  -11.2  eV  in  tie  emission 
spectra,  is  considerably  broader  than  either  of  the  nonbonding  bands 
because  of  the  larger  overlap  of  bonding  orbitals.  The  conduction 
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band,  which  begins  at  9.0  eV  above  the  top  of  the  valence  band, 

includes  the  Si-0  antibonding  orbitals. 

Figure  2  summarizes  schematically  the  electronic  structure  of  Si02 
and  hence  the  expected  optical  transitions  from  the  two  nonbonding 
orbitals  in  the  VB  to  the  unfilled  3d  level  of  Si,  which  can  be  probed 
by  no  reflectance  measurements.  For  bond  angles  less  than  180°,  due 
to  crystal  field  effects  the  degeneracy  of  p-orbitals  will  be  removed 
and  the  bottom  nonbonding  orbital  will  lower  its  energy  by  s  -  p 
orbital  mixing  (Fig.  2a).  The  VB  width,  AE,  shown  in  Fig.  2b  due 
to  crystal  field  splitting  alone  will  have  a  functional  dependence 
of  0  as  shown  in  Fig.  2c.  Figure  2  summarizes  the  expected  optical 
and  dielectric  properties  of  Si02  resulting  from  the  electronic  trans¬ 
itions  from  the  two  top  nonbonding  orbitals  of  the  VB  mainly  to  he 
unfilled  Si  3d-level  in  the  CB  as  a  function  of  0.  The  transitions 
mainly  control  AE  if  it  came  only  due  to  crystal  field  effects.  The 
uv  reflectance  spectra  should  probe  all  such  transitions  and  their 

dependence  on  0. 

Coesite  was  not  studied  due  to  the  unavailability  of  this  part¬ 
icular  allotrope  of  Si02  with  sufficient  macroscopic  size  and  optically 
good  surface  for  reflectivity  studies.  Synthetic  cristobalite  samples 
were  preposed  by  soaking  amorphous  quartz  discs  at  1400°C  for  90  hours 
followed  by  a  second  soaking  at  865°C  for  30  hours  before  furnace 
quenching  to  room  temperature.  On  the  basis  of  the  twelve  strongest 
reflections  found  by  x-ray  diffraction,  samples  were  identified  as 
100  percent  a-cristobalite  with  a  <101>  texture.  Tridymite  was  not 
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independently  studied  as  the  bond  angle,  6,  for  this  allotrope  is  the 

831,16  “S  ““  °f  -i-oballte.  The  fused  quartz  (optical  grade)  and 
the  Z-cut  a -quartz  samples  were  used  to  represent  amorphous  and 
crystalline  allotropes,  respectively.  For  reflectivity  measurements 
a  McPhe arson-2 25  type  monochrometer  was  used  which  Is  capable  of 
generating  uv  light  In  the  energy  range  0-15  eV.  Light  from  a 
thorlatad  tungsten  filament  using  a  glow  discharge  in  hydrogen  gas  was 
used  as  the  source  of  light.  The  light  goes  through  concave  gratings 
before  coming  to  the  exit  slit.  Optics  from  the  exit  slit  to  the 
sample  holder  assembly,  which  Is  vacuum  sealed  to  the  monochrometer, 
is  shown  in  Fig.  3  for  the  reflected  light  configuration.  The  incid¬ 
ent  light  intensity  configuration  was  attained  by  rotating  the  mirror 

about  the  axis  perpendicular  to  the  plane  of  the  paper.  Both  the 

incident  and  the  reflected  Intensity  were  detected  by  the  sodium 

salicylate  window  of  the  photomultiplier  tube.  In  these  particular 

optical  arrangements,  in  either  the  reflected  or  incident  light 

configurations,  the  pencil  of  light  hits  about  the  same  location 

of  the  mirror  only  once,  thus  minimizing  errors  in  the  reflectivity 
data. 

Figure  4  shows  the  reflectivity  data  taken  at  0.05  eV  intervals. 

The  amorphous  SiO.,  shows  two  distinct  peaks  located  at  10.3  and  11.55 
eV,  respectively,  corresponding  to  the  electronic  transitions  from 
the  two  nonbonding  orbitals  of  the  VB  edge  mainly  to  the  unfilled  Si 
d-level  in  the  CB.  The  width  of  the  VB  edge  4E,  is  found  to  be  1.25  eV. 
The  reflectance  spectra  for  the  crystalline  quartz  was  strikingly 
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similar,  as  observed  by  other a  previously,  excepL  that  -he  first  peak 
occurring  at  10.25  eV  is  much  sharper  and  that  in  the  second  peak, 
occurring  at  10.55  eV,  there  are  considerable  fine  structures.  This 
is  possibly  due  to  other  inter-band  and/or  intraband  transitions  frcm 
the  VB  edge  to  other  allowed  transitions  to  the  bottom  of  the  CB.  The 
results  for  crystobalite  are  somewhat  unusual  in  that  corresponding 
peaks  are  collapsed  around  10.7  eV  to  a  twin  peak  occurring  at  10.35 
and  10.9  eV,  respectively,  with  additional  fine  structures. 

In  spite  of  the  drastic  change  in  0,  reflectance  spectra  of  the 
three  allotropes  are  relatively  insensitive  to  0.  In  addition,  if 
the  width  of  the  VB  edge  were  determined  by  the  crystal  field  effects 
alone,  AE  would  have  collapsed  to  zero  for  cristobalite  and  a  single 
peak  would  be  expected,  which  was  not  the  case.  This  indicates  that 
AE  is  not  as  narrow  as  we  previously  thought  and  that  the  band  effects 
due  to  the  overlap  of  next -nearest  neighbor  oxygen-oxygen  orbitals 
is  rather  high  in  the  crystobalite  lattice. 

Kramer-Kronig  analyses  were  made  on  the  reflectivity  data  to 
find  the  functional  dependence  of  dielectric  (e^  e2)  and  other 
derived  optical  (n,  k)  properties  on  0.  and  e2  are  the  real  and 
imaginary  part  of  the  frequency-dependent  dielectric  function,  e(w) 
and  n  and  k  are,  respectively,  the  index  of  refraction  and  extinction 
coefficient,  (e  -  e±  +  ie.,;  (n  +  ik)?  -  ^  +  i^)  Figures  5  and  6 
show  the  variations  of  (e^  e.,)  and  (n,  k),  respectively,  as  a  function 
of  photon  energy  for  all  the  three  stated  allotropes.  They  show 
similar  variations,  as  seen  in  the  raw  reflectivity  data,  except  that 
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the  fine  structures  are  amplified.  The  variations  of  e ^  as  a  function 
of  photon  energy  from  Fig.  5,  below  the  absorption  edge,  can  be  best 
fitted  to  a  parabola  and  the  results  are  shown  in  Fig.  7.  The 
functional  dependence  of  /ej  near  the  VB  edge  on  the  photon  energy 
locates  exactly  the  VB  edge,  which  was  found  tp  be  9.05.9.1  and  9.25  eV 
for  cristobalite ,  amorphous  and  crystalline  quartz,  respectively. 

The  reflectance  spectra  are  relatively  similar  for  cristobalite, 
amorphous  and  crystalline  quartz.  This  leads  us  to  believe  a  simple 
model  based  on  the  oscillations  from  an  unhybridized  0  p  orbitals 
to  the  unfilled  Si  3d-level  in  the  CB. 

The  width  of  the  VB  edge  (AE)  is  comparable  for  all  three 
allotropes  and  hence  the  band  broadening  effects  due  to  the  overlap 
of  two  next-nearest-neighbor  oxygen  orbitals  is  high  enough  to  give 
a  finite  width  for  cristobalite  (0  «  180°),  where  there  is  no  crystal 
field  effects. 

AE  is  of  the  order  of  1  eV  which  is  broader  than  previously 
thought;  however,  it  is  narrow  enough  to  cause  localization  and 
trapping  of  holes  in  an  SiOj  lattice,  relatively  defect  free,  when 
SiO^  is  subject  to  an  ionizing  radiation. 
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—8 


I2i 


0° 


Figure  2:  Electronic  structure  and  the  valence  ban^  (VB)  edge  in  SiC^ 
(a)  For  0  <  180°,  the  mixing  of  lower  nonbonding  2p  orbital  with  2s 
orbital  due  to  crystal  field  effects  (b)  the  top  of  the  VB  edge  and  its 
width,  AE  (c)  the  functional  dependence  of  location  of  VB  edge  and  AE 
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B.  STRUCTURAL  TRANSFORMATIONS  AS  OBSERVED  BY  TEM  DURING 
ELECTRICAL  SWITCHING  IN  AMORPHOUS  Ge-Te 

R.  B.  Laibowitz,  S.  R.  Herd  and  P.  Chaudhari 

ABSTRACT:  Structural  transformations  were  induced  In  thin  films  of 
amorphous  Ge-Te  by  an  applied  voltage  while  the  sample  was  simultaneously 
viewed  via  transmission  electron  microscopy.  Three  resistance  states 
were  observed,  the  highest  for  the  completely  amorphous  phase,  and  the 
lower  two  for  states  consisting  of  crystalline  Te  in  both  an  amorphous 
material  and  in  crystalline  GeTe.  Switching  between  the  lower  two 
states  can  be  accomplished,  but  with  increasing  material  degradation. 

In  situ  electron  diffraction  was  used  to  identify  the  various  phases. 
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When  thin  films  of  amorphous  Ge-Te  are  sandwiched  between  conducting 
electrodes  and  biased,  electrical  switching  between  two  resistance  states 
occurs  (Ovshinsky  1968).  In  general,  two  distinct  bistable  resistance 
characteristics  have  been  observed:  one  in  which  the  resistance  state 
can  be  maintained  at  zero  bias  (memory  state)  and  a  second  in  which  a 
holding  voltage  is  necessary  to  maintain  the  particular  resistance. 

Whi.e  work  has  been  done  attempting  to  understand  these  properties  in 
Ge-Te  and  other  amorphous  materials  (Bosnell  and  Thomas  1973),  the 
physical  mechanisms  responsible  for  the  conductivity  changes  remain  as 
yet  unresolved.  This  is  especially  true  of  the  morphological  trans¬ 
formations  that  occur  in  these  systems.  In  this  letter,  we  report  on 
the  results  of  an  investigation  in  which  we  determined  the  microscopic 
structural  changes  in  the  amorphous  Ge^Te^  during  switching  and 
related  these  changes  to  the  presence  and  behavior  of  the  memory  state. 

The  technique  we  used  was  to  directly  view  tha  amorphous  material  via 
high  resolution  transmission  electron  microscopy  (TEM)  while  simultan¬ 
eously  taking  I-V  curves  which  showed  the  switching  behavior  (Chaudhari 
and  Laibowitz  1972;  Laibowitz,  Chaudhari,  Herd  and  Mader  1973).  In 
addition,  in  situ  electron  diffraction  was  used  to  identify  the 
structure  after  phase  changes  ware  observed. 

The  experimental  set-up  will  be  briefly  reviewed  here;  a  more 
complete  discussion  of  the  sample  configuration  can  be  found  in  our 
earlier  report  (Chaudhari  and  Laibowitz  1972).  Basically,  two  closely 
spaced  metal  electrodes  were  fabricated  on  the  surface  of  an  oxidized 
Si  wafer  by  electron  beam  techniques  (Hatzakis  and  Broers  1971).  The 
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width  of  the  gap  between  the  two  electrodes  could  be  varied  from  about 
O.ly  to  several  microns.  Amorphous  Ge-Te  of  approximately  Ge1(.Teoc 
composition  was  then  coevaporated  from  separate  sources  onto  this 
sample,  filling  the  gap.  The  substrate  beneath  the  gap  was  then  sel¬ 
ectively  back-etched  until  the  remaining  SiC^  with  the  Ge-Te  on 
top  was  transparent  to  100  KeV  electrons  used  in  the  microscope.  In 
general,  the  thickness  of  the  Ge-Te  was  about  O.ly  or  less  with  a 
residual  Si02  thickness  of  about  0.02y.  Leads  were  then  connected  to 
the  metal  electrodes  (Al)  and  the  sample  mounted  in  the  microscope. 

Either  ramp  or  pulsed  biases  could  be  applied  to  the  sample  by  external 
sources  and  the  response  of  the  sample  was  monitored  on  an  oscilloscope. 

When  a  bias  was  first  applied  across  the  gap,  i.e.,  across  the 
Ge-Te,  a  very  high  resistance  state  was  obtained.  Upon  increasing  the 
voltage,  a  partial  breakdown  (or  forming)  occurred  resulting  in  a  low 
resistance  state.  This  low  resistance  state  is  caused  by  conduction 
through  crystalline  filaments  (Chaudhari  and  Laibowitz  1972)  which 
bridge  the  gap  as  shown  in  Fig.  1.  However,  this  type  filament  could 
only  be  observed  in  thicker  (>  O.ly)  Ge-Te  films  which  allowed  visual 
observation  but  were  too  thick  for  electron  diffraction. 

Subsequent  measurements  on  thinner  (0.06y)  Ge-Te  films  allowed 
analysis  by  electron  diffraction  and  are  shown  in  Figs.  2a  and  2b.  Again, 
after  the  low  resistance  state  was  obtained,  we  found  crystallization, 
although  the  morphology  of  the  crystals  in  the  thin  films  is  somewhat 
different.  The  speckled  regions  of  Figs.  2a  and  2b  show  crystalline  Te 
particles  in  an  amorphous  Ge-Te  matrix  which  are'  similar  to  the  crystals 
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observed  after  the  first  forming  took  place.  Upon  applying  a  high 
current  pulse,  the  Te  redisaolves  and  the  material  becomes  amorphqus 
again,  i.e.,  the  material  attempts  to  return  to  its  initial  high 
resistivity  state.  Some  damage  usually  occurs  in  these  films  during 
switching  with  a  high  current  pulse  due  to  material  flow. 

Figure  2a  shows  the  sample  in  its  lowest  resistance  state  after  a 
few  switching  cycles.  We  find  a  band  of  small  Te  crystals  as  observed 
after  the  initial  forming.  In  addition,  near  the  top  of  Fig.  2a  where 
the  material  thickened  due  to  flow  during  the  previous  resetting,  we 
find  two  kinds  of  crystallization.  Conduction  occurs  through  both  Te 
dendrites  in  amorphous  Ge-Te  (right  marker)  and  Te  in  crystalline  GeTe 
(left  marker).  Figure  3  shows  the  selected  area  diffraction  pattern 
and  micrographs  of  this  region  at  a  higher  magnification.  Most  of  the 
spotty  rings  in  the  selected  area  diffraction  pattern  (Fig.  3a)  are  due 
to  randomly  oriented  Te  crystals  (hex.,  aQ  -  4.46  X;  cq  -  5.93  X).  The 
short  arced  spot  marked  by  a  circle  in  the  diffraction  pattern  is  due 
to  the  <200>  reflection  of  GeTe.  A  dark  field  micrograph  (Fig.  3b) 
taken  with  this  reflection  shows  that  the  Te  crystals  visible  in  bright 
field  in  Fig.  3c  are  embedded  in  essentially  a  singla  crystal  of  GeTe. 
There  are  indications  that  the  GeTe  is  in  the  rhorabohedral  form  (a  - 
5,98  X)  (a  -  88.35°),  (Chopra  and  Bahl  1969)  however,  positive  distin¬ 
ction  from  the  cubic  phase  (a  01  5.86  X)  is  difficult  in  this  case. 

Subsequent  switching  with  a  high  current  pulse  rendered  this  same 
area  amorphous  again  as  shown  in  Fig.  2b  and  simultaneously  the  sample 
returned  to  a  higher  resistance  state.  Temperatures  reached  during 
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switching  are  thus  high  enough  to  melt  both  crystal  phases  and  quench 
the  melt  back  into  an  amorphous  solid.  Again,  material  flow  took  place 
during  this  resetting  operation  indicated  by  the  thick  ridge  of  amorphous 
material  in  the  top  section  of  Fig.  2b.  Also  visible  is  a  sharp  demark- 
ation  line  where  the  temperature  was  not  sufficient  to  redissolve  "he 
small  Te  crystals  in  the  lower  half  of  the  micrograph. 

Thus,  our  experiment  which  combines  transmission  electron  microscopy 
with  electrical  switching  has  shown  that  three  conductivity  states  can 
be  produced  in  Ge1cTeoc  and  that  these  three  states  have  specific 
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morphologies  associated  with  them.  The  completely  amorphous  phase  has 
the  highest  resistance  while  the  combination  of  Te  crystallites  in  an 
amorphous  background  has  a  lower  resistance.  The  lowest  resistance 
state  is  achieved  with  both  crystalline  GeTe  and  crystalline  Te  in 
the  conducting  filamentary  region.  While  it  is  possible  to  cycle  (or 
switch  several  times  between  these  resistance  states,  the  damage  in 
these  very  thin  films  due  to  material  flow  usually  increases  leading 


to  eventual  bum  out. 
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Crystalline  filament  formed  in  thick  (>0.1p)  amorphous 
Ge-Te  matrix  30.000  x  Mag. 
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Figure  3a: 


Selected  area  diffraction  of  crystalline  Te  and  GeTe. 
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Figure  3b: 


Dark  field  mircograph  of  GeTe  crystal  taken  with  <200> 
reflection  as  indicated  in  Fig.  3a  Mag:  200, 000  x 
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Ill,  INSTABILITY  ASSOCIATED  WITH  METAL-GLASS  INTERACTIONS 
In  the  first  semi-annual  report,  we  have  reported  a  study  of  the 
interaction  between  vanadium  and  thermally-grown  SiO,,.  The  interaction 
decomposes  Si02  and  results  in  the  formation  of  a  vanadium-rich 
silicide,  V3Si.  It  was  suggested  in  the  report  that  the  metal-glass 
interaction  may  not  be  unique  to  V,  but  is  a  rather  general  phenomenon 
between  glass  and  transition  metals  that  form  highly  stable  oxides. 

The  transition  metal  oxide  which  is  more  stable  than  SiO,,  provides 
he  driving  force  for  the  reactions.  In  a  continuous  effort  to  under¬ 
standing  the  glass-metal  interaction,  we  have  further  investigated  the 
Interaction  between  oxidized  silicon  and  Ti  and  Nb  using  a  combined 
technique  of  x-ray  diffraction  and  ion  backscattering.  Indeed  we 
found  that  the  results  agree  very  well  with  our  expectation  that  both 
Ti  and  Nb  decompose  Si02  at  temperatures  above  700°C  and  900°C, 
respectively  and  form  metal-rich  silicides,  Ti5Si3  and  Nb^i.  We  note 
that  the  Nb3Si  obtained  is  not  a  high  Tc  superconducting  phase;  it 
has  a  F.C.C.  structure  the  same  as  the  ordered  Cu3Au. 

Silicide  formation  in  silicon  technology  is  most  commonly  obtained 
by  a  direct  reaction  between  silicon  and  a  vapor  deposited  metallic 
film.  To  compare  our  glass-metal  reaction  to  silicon-mt ;al  reaction 
from  the  point  of  view  of  silicide  formation,  we  have  studied  contact 
reaction  between  Ti  and  Si  and  also  Nb  and  Si.  The  results  again 
agree  with  our  earlier  findings  in  the  case  of  V  and  Si;  the  silicide 
formed  by  metal-silicon  reaction  is  silicon-rich  TiSi2  and  NbSi2> 
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There  exists  a  distinct  difference  between  metal-silicon  and 
metal-glass  reactions  from  the  point  of  view  of  silicide  formation. 

As  a  first  step  to  understanding  the  cause  of  the  difference,  and 
also  because  there  has  recently  appeared  a  large  quantity  of  literature 
on  silicide  formation,  we  have  made  a  review  of  silicide  formation 
from  the  viewpoint  of  their  structure,  kinetics,  and  stability.  The 
review  is  intended  to  update  our  understanding  of  this  subject  and 
also  to  suggest  future  research  directions,  in  particular,  we  have 
pointed  out  that  marker  motion  in  silicides  and  the  effect  of  a  thin 
oxide  on  metal-silicon  interactions  are  both  crucial  areas.  Also 
included  in  the  review  is  a  critical  evaluation  of  sample  preparation 
and  experimental  techniques  of  ion  backscattering  and  glancing  angle 
x-ray  diffraction  as  they  are  applied  to  study  thin  film  interactions. 
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A.  ANALYSIS  OF  THIN-FILM  STRUCTURES  WITH  NUCLEAR  BACKSCATTERING 
AND  X-RAY  DIFFRACTION 

J.  W  Mayer  and  K.  N.  Tu 

ABSTRACT:  Backscattering  of  MeV  He  Ions  and  Seemann-Bohlin  x— ray 
^ raction  techniques  have  been  used  to  study  sillclde  formation  on 
Si  and  SiC^  covered  with  evaporated  metal  films.  Backscattering 
techniques  provide  information  on  the  composition  of  thin  film 
structures  as  a  function  of  depth.  The  glancing  angle  x-ray  technique 
provides  identification  of  phases  and  structural  information.  Examples 
are  given  of  V  on  Si  and  on  SiC^  to  illustrate  the  major  features  of 
these  analysis  techniques.  We  also  give  a  general  review  of  recent 
studies  of  silicide  formation. 


I  INTRODUCTION 


The  formation  of  silicide  layers  has  played  an  important  role  in 
integrated  circuit  technology.  Ohmic  contacts  and  barriers  on  Si  are 
formed  by  evaporating  a  metal  layer  (Pd,  Pt,  Ti,  ...)  on  Si  and  subsequently 
heating  to  form  the  metal-silicide . ^ ^  We  have  utilized  backscattering 
techniques  with  MeV  He  ions  to  study  the  kinetics  of  the  silicide  growth 
and  glancing  angle  x-ray  diffraction  techniques  to  identify  the  phases 
that  are  formed.  This  paper  is  a  review  of  the  capabilities  of  these 
analytical  techniques  and  of  the  results  obtained  in  our  and  other  labora¬ 
tories  concerning  metal-silicon  interactions. 

Hie  process  steps  of  evaporation  and  heat  treatment  are  simple  and 

easily  adapted  to  device  production.  The  factors  which  govern  the  silicide 

formation  are  more  complex.  Of  the  twelve  metal/silicon  systems  studied 
(2-17) 

in  some  detail  ,  all  have  between  three  and  seven  silicide  structures 

(18—21) 

identified  in  bulk  samples  .  In  thin  film  structures,  typically  only 

one  and  in  a  few  cases  two  silicide  phases  are  formed.  The  lowest  eutectic 
temperatures  given  in  the  phase  diagrams  range  between  740°C  and  1400°C, 
yet  the  silicides  form  in  thin  film  samples  at  temperatures  between  200° 
and  800°c .  In  most  cases  the  formation  occurs  at  temperatures  of  about 
one-third  to  one-half  the  melting  point  (in  °K) .  This  is  a  sure  sign  that 
solid-solid  rather  than  solid-liquidus  reactions  are  involved.  It  has 
recently  been  found  that  silicide  layers  are  formed  during  heat  treatment 
of  metal  layers  on  Usually  a  different  phase  is  formed  and 

the  reaction  temperature  is  100  to  200 °C  higher  than  for  the  metal/silicon 
sys  terns . 
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Prom  an  operational  standpoint  it  is  of  interest  to  determine  the 
temperatures  at  which  silicide  growth  occurs,  growth  rates,  identity  of 
the  phases,  their  stability,  and  the  influence  of  the  interface  conditions 
(presence  of  oxide  layer) .  These  aspects  are  of  concern  in  integrated 
circuit  technology  and  can  be  answered  directly  by  the  combined  use  of 
MeV  He  backs cat tering  and  glancing-angle  x-ray  diffraction  techniques. 

The  underlying  problem  is  to  determine  what  factors  govern  the  nucleation 
and  growth  of  the  silicide.  Some  insight  can  be  gained  by  measuring 
the  kinetics  of  the  process  and  determining  which  species,  Si  or  metal, 

diffuses  through  the  silicide  layer.  Again,  backscattering  techniques 
can  be  used. 

II  ANALYTICAL  TECHNIQUES 
4 

A.  MeV  He  Ion  Backscattering 

This  technique  has  been  described  in  detail  previously ^ ~2~2^  and 
will  be  discussed  only  briefly.  As  shown  in  Fig.  1,  the  sample  is  mounted 
in  a  chairher  (modest  vacuum  requirements,  £  10~5  -  lo"6  Torr)  and  bombarded 
with  monoenergetic  He  ions  (typically  at  2  MeV)  at  current  levels  of  20-50  na 
over  about  a  2  mm  x  2  rrm  area.  Only  a  small  fraction  of  the  incident  particles 
are  backscattered  so  the  beam  is  to  first  order  not  attenuated.  The  silicon 
surface  barrier  detector  produces  a  voltage  pulse  that  is  proportional  to 
the  energy  of  the  individual  backscattered  He  ions.  The  pulses  are  amplified 
and  stored  in  a  pulse-height  analyzter  and  at  the  end  of  the  analysis 

(generally  15-30  min)  the  energy  spectrum  of  backscattered  particles  is 
displayed. 
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The  energy  of  particles  scattered  from  surface  atoms  is  determined  by 

the  mass  of  target  atom  through  scattering  kinematics  K  (billiard-ball 

kinetics) ,  the  numer  of  particles  backscattered  by  the  scattering  cross- 

section  o  (Rutherford  scattering  with  o  a  (Z/E) 2  in  this  energy  region) 

the  energy  loss  in  traversing  the  film  by  stopping  cross  section  e. 

these  parameters,  K,  o  and  £,  are  well  known  with  the  major  uncertainty 

in  the  stopping  cross-section  e  where  the  values  are  generally  established 
(24) 

within  +  5%  .  Consequently,  these  backscattering  spectra  give  directly 

and  quantitatively  the  distribution  of  elements  as  a  function  of  depth  in 

the  target.  In  effect  it  provides  mass-sensitive  depth  microscopy. 

It  should  be  noted  that  the  depth  scale  is  determined  by  the  number 
2 

of  target  atoms/cm  in  an  incremental  layer  not  by  the  thickness  of  the 
layer.  Cast  in  other  terms,  the  density  of  a  thin  film  can  be  determined 
from  measurement  of  the  number  of  target  atoms/cm2  (or  from  nhe  difference 
in  energy  between  particles  scattered  from  the  front  and  rear  surface  of  a 
film)  and  from  an  independent  measurement  of  the  film  thickness  by  inter- 
f e rome trie  techniques.  For  simplicity,  the  film  thicknesses  are  usually 

e 

given  in  dimensions  of  A  with  the  assumption  of  bulk  density. 

The  advantages  of  backscattering  then  are:  1)  quantitative  measurement 
of  composition,  2)  non-destructive  determination  of  depth  profiles,  and 
3)  simple  and  direct  interpretation  of  experimental  data.  The  measurements 
themselves  are  quite  easy  to  preform. 

Some  of  these  features  are  shown  in  Fig.  2  which  gives  backscattering 
spectra  for  an  evaporated  layer  of  V  on  Si  before  and  after  heat  treatment. 
For  the  as-deposited  case  (dashed  line)  ,  the  broad  peak  is  the  component 
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of  spectrum  due  to  particles  backscattered  from  V.  The  measured  energy 
difference  between  particles  scattered  from  the  front  surface  and  the  back 
surface  of  the  V  layer,  i.e.  the  width  of  the  peak,  is  300  keV  which 
corresponds  to  2.24  x  1018  V  atoms/cm2  or  a  thickness  of  3100A  if  bulk 
density  is  assumed.  After  heat  treatment  a  step  appears  in  both  the  V 
and  Si  portions  of  the  spectrum.  This  step  corresponds  to  the  formation 
of  a  f-ilicide  layer.  The  composition  can  be  determined  from  the  ratio  of 
the  heights,  H,  of  the  V  and  Si  steps,  here  1.35,  when  corrected  for  the 
ratio  of  scattering  cross  sections,  ^si/ov  -  0.36.  This  gives  a  V/Si  ratio 
of  .485  which  must  be  corrected  upwards  by  a  factor  of  1.04  because  of  the 
differences  in  stopping  cross-section  for  particles  scattered  from  V  and 
from  Si.  This  result  shows  that  in  the  layer  there  are  two  Si  atoms  for 
every  V  atom,  an  indication  that  the  compound  is  formed.  However, 

positive  identification  requires  structural  data  which  is  provided  by 
x-ray  diffraction  (see  B,  below). 

The  thickness  t  of  the  composite  layer  is  determined  from  the  energy 

width  Ae,  here  AE  “  140  keV.  Assuming  bulk  density  this  gives  a  thickness 
0 

of  1700A.  The  growth  kinetics  of  the  composite  layer  can  be  found  from 
the  time  dependence  of  Ae  for  various  heat  treatment  cycles. 

The  fact  that  the  edges  of  the  steps  are  sharp  in  Fig.  2  demonstrates 
that  the  layer  is  uniform  in  thickness  over  the  dimensions  of  the  beam. 
Laterally  non-uniformities  introduce  a  broadening  in  the  edges  and  hence 
uncertainties  in  thickness  determination.  If  local  cracking  or  flaking 
of  the  film  occurs  during  thermal  processing,  backscattering  spectra  are 
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misleading.  Hence,  one  should  also  gamine  lateral  »ifor»ity  by 

optical  micros  copy  or  preferably  by  use  of  a  scanning  electron  microscope. 

Backs catte ring  spectra  are  simple  to  interpret  when  there  are  o  y 
few  elements  in  the  target.  For  multilayer  targets,  the  spectra 

complex .  Figure  3  shows  spectra  for  nortal  incidence  ana  45*  beam,  to 

,  _  mur  _7nc  reflective  coating.  The 

target  orientation  for  a  twelve  layer  ThF4 

.  lavers  and  the  valleys  show 

distinct  peaks  reflect  the  Th  in  the  twelve  ThF4  Y 

me  position  and  thinness  of  the  ZnS  layers.  TOO  contribution  the 

s,  and  F  can  not  directly  be  seen  in  the  spectra.  The  arr»s  indicate  their 
energy  position  it  located  at  the  snrFacc.  The  too  pmaK.  at  the  highest 
energy  are  flat  topped  and  have  sharp  edges  indicated  well  defined  lay 
of  uniform  coalition.  *  l—  energies  (greater  depths  in  - 
structure)  the  broadening  OF  the  edges  reflect  both  a  loss  of  depth  resolntlon 
due  to  energy  straggling  and  interference  from  lower  mass  elements. 


B.  Seemann-Bohlin  X-ray  Diffraction 

The  Seemann-Bohlin  configuration  for  structural  analysis  of  thin 

ffraetion  using  a  focused  beam 
films  is  basically  a  glancing  angle  x-ray  diffr 


uith  a  fixed  angle  of  incidence,  it  is  hno«,  that  to  study  the  stmeture 


Of  a  polycrystalline  thin  film  of  thickness  of  the  order  of  10  to  1 

•  _ i-  MttoVi  matAfl 


UL  a  J -  . 

difficult  by  conventional  x-ray  techniques  because  there  is  not  much  materia 
fur  scattering.  Hoveve,  scattered  beams  can  be  strengthened  by  increasing 
the  path  of  ray  in  the  film  by  *e  use  of  a  fixed  iuerdent  beam  at  a  near- 

i  c  \  yaA  incident  be  sin  c&n  be 
grazing  incidence,  mi  »e  note  that  such  a  taxed  rnerd 

v.  *  ran  fhB  focusing  Seemann-Bohlin  configuration 
obtained  in  a  diffractometer  based  on  the  focusing 

I  ZD)  ^/inaina 


uuuaAncv*  —  —  /o  c\  _  .  _ 

'2t>)  In  the  focusing 

but  not  on  the  conventional  Bragg-Brentano  configuratron 
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Seeroann-Bohlin  configuration,  see  Fig.  4,  the  specimen  is  placed  on  the 
circumference  of  the  diffraction  circle  and  the  angle  of  incidence  is  fixed 
and  can  be  made  as  small  as  a  few  degrees,  for  example,  an  incident  angle  of 
6.4°  in  this  case  will  increase  the  length  of  the  path  of  the  beam  in  the 
specimen  to  about  9  times  its  thickness.  The  specimen,  the  focus  "F"  of  the 
monochromatic  incident  beam  and  the  focus  of  the  diffracted  beams  all  lie 
ideally  on  the  circumference  of  the  diffraction  circle.  Based  on  this  principle, 

/  2f>\ 

Feder  and  Berry  have  designed  and  built  a  x-ray  diffractometer  with  a 
diffraction  circle  of  20  inches  in  diameter.  An  important  feature  of  the 
diffractometer  is  that  it  employs  a  pyrolytic  graphite  monochromator  crystal 
to  produce  a  high  intensity  monochromatic  Cuka  radiation  as  the  incident 
beam.  The  intensity  obtained  from  the  graphite  crystal  is  15  times  greater 
than  that  obtained  from  a  LiF  crystal.  Also  a  circular  helium  chamber  of 
18  inches  in  diameter  is  interposed  in  the  path  of  the  diffracted  beams  to 
reduce  the  atmospheric  absorption. 

(26) 

The  performance  of  the  diffractometer  has  been  discussed  elsewhere 
VJe  mention  briefly  here  that  it  has  detected  a  polycrystalline  copper  film 

O 

of  150A  and  shewed  a  precision  lattice  parameter  measurement  of  a  nunber  of 

o  0 

1000A  nickel  films  with  an  accuracy  of  +_  0.0001A,  i.e.  1  part  in  35000.  In 

addition  to  the  high  sensitivity  of  detecting  phases ,  the  Seemann-Bohlin 

diffraction  allows  convenient  structural  analysis  of  grain  size,  biaxial 

strain  and  the  distribution  of  pole  densities.  This  is  because  its  reciprocal 

lattice  vector  "H"  of  Bragg  d  ffractioris  are  not  parallel  to  the  normal  "n" 

of  the  film  surface  but  rather  make  inclination  angles  with  the  normal. 

The  angle  increases  with  the  order  of  the  reflections.  Thus  each  reflection 
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measures  in  a  different  direction  the  lattice  parameters,  the  size  and 
the  total  amount  of  coherently  diffracted  regions.  It  is  clear  that  this 
type  of  diffractometer  can  deliver  a  great  deal  of  structural  information 
of  polycrystalline  thin  films.  We  should  note  that  single  crystal  filiu^ 
are  not  suitable  for  Seemann-Bohlin  diffractometer  and  the  diffraction 
patterns  may  show  no  peaks  at  all.  This  may  be  regarded  as  advantageous 
for  a  polycrystalline  film  on  a  single  crystal  substrate  such  as  a  metal 
film  on  a  Si  wafer  since  we  do  not  need  to  worry  about  the  interference 
of  Si  reflections. 

It  seems  that  for  a  structural  study  of  thin  films  of  thickness  around 

O 

1000A  and  less,  we  can  use  transmission  electron  microscope,  for  films 
thicker  than  several  microns  we  can  use  conventional  x-ray  equipment,  and 
for  films  in  between,  we  use  Seemann-Bohlin  type  x— ray  diffractometer. 

The  diffractometer  is  expected  to  be  very  powerful  in  the  investigation 

(27) 

of  the  interaction  in  thin  films  ,  because  the  reaction  products  can  be 
identified  by  their  x-ray  reflections  very  early  in  the  diffusion  cycle. 

O 

Figure  4  shews  the  x-ray  pattern  of  a  sample  of  3200A  of  vanadium 
on  a  Si  wafer  that  has  been  heat  treated  at  800°C  for  15  minutes.  The 
pattern  was  obtained  by  scanning  the  sample  with  CuKa  radiation  at  steps 
of  0.15*  (40)  increment  and  with  a  counting  time  of  30  seconds  at  each 
step.  The  peaks  have  been  indexed  as  reflections  of  VS^  according  to 
AS7M  card  #13-260.  In  this  sample,  the  x-ray  diffraction  data  confirm 
the  composition  analysis  determined  from  backscattering  data  (Fig.  2) . 
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Ill  ANALYSIS  OF  SILICIDE  FORMATION 
A.  General  Precautions 

One  of  the  objectives  of  this  tvDe  ,  • 

growth  kinetics  it  is  ^  ^  *° 

ia  th  d  .  “  *  preparation 

is  the  dominant  factor.  ^ere  arP  *.K 

jere  are  three  areas  which  mnqt  v  Q 

wmcn  must  La  considered: 

interface,  contamination,  and  stress. 

/'  i£SS£Ji££'  tte“i"^‘«„„llllOTlirflcBititb( 
considered  when  preparing  samples  it  k 

sl  ...  •  t  has  been  shown  that  thin  oxide  layer 

r  *  ret*rd  “  .  -  -i.  jl 

Z  —  —  orientations  ^28)  ,„d  _ 

*  0rlentSti0n  Playa  a  r°le-  *»«“«  *»  HP  ana  then  rinsing  prior  to 
™9  in  the  evaporation  does  not  „ee  „  ^  ^ 

-race,  one  approach  is  to  sputter  ciean  the  surface  a„d  then  aeposit 

SIT  T'  "ith°Ut  brelk1"9  ^  “*-*  ^  ~  reproducible 

surface  conditions,  the  sinPnn 

SUrf3Ce  is  he^ily  damaged  and  the 

SDllffavta. _ . 


sputterin,  gas  is  retained  in  th.  surface  layer. 


(29) 


*•  Co^Hinatio,.  Oxygen  and  other  anient  .intonations  can  be 

“  rcd"“d  d“ring  “  Processing.  Trace  inclusion  of 

oxygen  are  hard  to  detect  h«  i, 

M  .  7  techniques ,  aithough  relatival 

massive  concentrations  (10  to  20  at 

20  at.%)  can  be  readily  identified.  There 

;  ^  bee"  3  P03Ui''e  -  -  o°nt  ami  nation  oith  retardation 

T:;hi  h°—  — -  -  ~  that  this  nay  be  a 

real  problem. ( 30 } 

3‘  —  °"e  °£  0,6  factors  to  overtop  is  fila 

rrachm,  or  peeling  due  to  stress  created  during  them,.!  processing. 
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Backscattenng  measurements  on  such  samples  are  unreliable.  Film  peeling 
is  often  enhanced  for  thicknesses  greater  than  a  few  thousand  Angstroms. 

In  some  cases  ,  this  problem  can  be  alleviated  by  maintaining  an  oxide  free 
interface  or  depositing  the  film  at  elevated  temperatures. 

There  is  of  course  the  polycrystalline  nature  of  the  deposited  and 
reac  ccl  film  that  should  be  considered.  Grain  size  and  grain  boundary 
diffusion  may  be  important  particularly  at  low  process  temperatures 
(T  <  400°C) .  These  problems  have  not  been^properly  addressed  in  silicide 
formation,  but  are  known  to  occur  in  metal-metal  thin  film  studies.  At 
higher  temperatures  bulk  diffusion  will  probably  dominate  over  grain 
boundary  effects. 


B.  An  Example:  V-Si  System 

Back.scattering  analysis  have  been  applied  to  nearly  a  dozen  silicide 
systems.  Two  of  which  (Pd-Si ( 3 '6 '9 '10 '12_14)  ^  Hf-si{7'8))  have  been 
studied  in  detail.  We  have  chosen  the  vanadium  silicide  system  as  an 
example  because  tie  problems  are  still  fresh  (as  the  study  has  not  been 
completed)  and  also  we  cat  consider  aspects  of  silicide  formation  on  SiC>2. 
Our  work  with  Ti,  Nd  and  V  deposited  on  Si02  shows  that  silicides 

are  formed  but  at  higher  processes  temperature  than  for  Si/metal  systems 
and  that  the  silicide  is  metal  rich. 


Figure  6  shows  backscattering  spectra  for  V  films  deposited  on  Si02<30). 
For  the  as-deposited  case  (dashed  line  in  Fig.  6B)  the  spectrum  shows  the 
contribution  from  V,  Si  in  the  oxide  layer  and  substrate,  and  oxygen  in 
the  oxide  layer.  The  oxygen  contribution  is  superin$>osed  on  that  from  the 


underlying  Si  substrate.  The  arrows  indicate  the  energy  position  of  these 
elements  when  at  the  surface.  Following  heat  treatment  in  oxygen  at  400°C 
there  is  a  change  in  the  composition  of  the  sample.  Oxygen  is  now  present 
at  the  outermost  layers  (I)  as  evident  in  Fig.  6A  from  the  step  in  the  V 
component  and  the  appearance  of  a  broad  mesa  whose  high  energy  edge  corresponds 
to  the  position  of  oxygen  at  the  surface.  There  is  not  a  significant  amount 
of  oxygen  in  underlying  V  layer  (II)  as  indicated  by  i)  the  height  of  the 
V  portion  of  the  spectrum  (layer  II)  coincides  with  the  as-deposited  case 
and  li)  there  is  a  dip  in  the  oxygen  component  of  the  spectrum  in  the  region 
between  that  from  oxygen  in  SiC>2  (III)  and  in  Vandium  oxide  (I)  . 

The  composition  of  the  Vandium  oxide  can  be  estimated  from  ratio  of 
the  heights  of  Vandium  and  oxygen  components  in  region  I  to  be  close  to 
that  of  V205.  Positive  identification  was  made  from  glancing  angle  x-ray 
diffraction  data  (Fig.  7).  The  peaks  correspond  to  those  indexed  for  V^. 

Referring  again  to  Fig.  6,  it  is  evident  that  there  is  a  dramatic 
change  in  composition  when  a  similar  sample  is  heated  to  800°c  in  vacuum. 

The  composition  of  the  layers  can  be  deduced  from  the  spectrum  to  3 how: 

Lajrer_I:  Approximately  equal  amounts  of  Vandium  and  oxygen  with  no 
appreciable  concentration  of  Si. 

— y~r  11 :  Three  times  as  much  V  as  Si  and  no  appreciable  oxygen. 

La^er  III:  SiC>2  but  of  reduced  thickness.  The  amount  of  oxygen 
removed  from  the  initial  thickness  of  Si02  (III)  corresponds  to  the  amount 
added  to  the  vanadium  oxide  layer  (I)  indicating  that  the  amount  of  oxygen 
in  the  composite  layers  is  conserved. 


The  x-ray  diffraction  data  for  this  sample  coincides  with  that  given 

by  Tu  et  al^1*  for  a  similar  sample  and  shows  the  presence  of  V3Si,  V5Si3 

and  V  0  .  The  major  feature  is  the  predominance  of  V3Si  for  this  sample 
2  5 

whereas  VSi2  is  formed  in  the  V/Si  samples. 

There  are  two  unusual  features  in  this  comparison  between  backscattering 

and  x-ray  diffraction  data.  Backs  cattaring  data  indicate  a  1:1  ratio  of  V 

to  0  and  diffraction  data  indicate  V^.  This  may  be  due  to  the  fact  that 

the  V  0  eutectic  is  at  a  lower  temperature (l8)  than  the  800°C  process 
2  5 

temperature.  We  suggest  that  oxygen  dissolves  in  the  vanadium  at  the  process 
temperature  and  then  the  V^,.  nucleates  as  crystallites  imbedded  in  a  matrix 
of  vanadium  when  the  sample  cools  below  the  eutectic.  However,  the  average 
V/o  ratio  over  the  entire  vanadium-oxide  film  thickness  could  be  near  unity 
as  indicated  by  backscattering  data.  The  other  question  concerns  the  location 
of  the  V5Si3  identified  in  x-ray  diffraction.  Backscattering  data(16)  do  not 
show  appreciable  amounts  of  Si  near  the  surface  (region  I  in  Fig.  6B)  and  the 
spectrum  height  ratio  for  the  silicide  layer  (region  II)  corresponds  to  V3Si. 
At  present  we  suggest  the  V,.Si3  is  located  near  the  interface  between  region 

I  and  II. 

This  example  clearly  shows  the  drastic  effect  that  a  thin  oxide  can 
have  on  both  growth  kinetics  and  silicide  exposition.  In  this  case,  the 
ratal-rich  silicide  is  of  interest  in  its  own  right  as  a  superconductor 

,  n.Jn) 

with  a  transition  temperature  of  17  K 
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IV  SILICIDE  FORMATION 


In  this  section  we  will  review  experimental  data  on  silicide  formation 
from  the  standpoint  of  structure,  thermodynamics  and  kinetics.  We  also 
include  some  speculations  and  suggestions  for  future  work.  Figure  8  gives 

a  general  impression  of  results  obtained  and  more  details  are  given  in 
Table  I. 


A.  Structure 

The  typical  silicide  which  is  formed  with  metal/Si  interaction  is  the 
Phase  which  is  most  rich  in  Si,  i.e.  PtSi  and  CrSi.,.  This  is  the  end  point 
which  one  expects  from  an  equilibrium  point  of  view.  There  should  be  two 
Phases,  Si  and  the  most  Si  rich  silicide.  It  is  surprising  that  more 
phases  do  not  appear  as  one  would  expect  all  phases  to  be  present  during 
the  reaction  stage.  Indeed  intermediate  phases  (Pd2Si,  Pt2Si,  HfSi  and 
NiSi)  have  been  detected  in  a  few  systems,  other  phases  may  be  present  but 
their  extent  is  determined  from  kinetic  considerations  of  diffusion  rate 
and  nucleation.  KidsonJ2  had  earlier  pointed  out  that  if  the  rate  constants 

of  diffusion  in  a  given  phase  are  small  the  phase  layer  may  be  too  thin  to 
be  detected. 

The  microstructure  of  HfSi,  HfSi,,,  VSi2,  PtSi  and  Pd^i  has  been 
measured  by  x-ray  diffraction.  As  a  first  approximation  the  grain  size  of 
the  silicide  is  about  the  same  ks  the  deposited  metallic  film  (£  200-500A) . 
There  have  not  been  serious  attempts  to  measure  strain  in  the  silicide 
layer  and  its  dependence  on  heat  Vestment  temperature.  This  is  clearly 
one  area  that  should  be  explored  further. 
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Even  though  the  crystal  structure  of  most  silicides  are  different 
from  Si,  they  tend  to  show  some  preferred  orientation.  The  most  striking 
case  is  Pd2Si  which  has  been  shown  by  both  reflection-electron  diffraction 
and  MeV  He  channeling  techniques  to  grow  epitaxially  on  <111>  si-10*12-14 
In  this  case  the  hexagonal  basal  plane  of  Pd^i  matches  quite  well  with 
the  <111>  Si  structure.  Similar  considerations  of  lattice  matching  suggest 
that  NiSi2,  CoSi2  and  FeSi,,  may  also  form  epitaxially.  This  conclusion 
has  not  been  tested  and  is  an  obvious  area  for  further  study. 

B.  Thermodynamics 

Silicide  formation  energies  have  been  tabulated33  for  the  following  four 
systems:  Mo/Si,  Ta/Si,  Ti/Si,  and  Zr/Si .  In  the  Mo/Si  system,  MoSi.,  is 
the  most  stable  phase  among  the  Mo-silicides  and  was  found  to  be  the 
predominant  phase  in  backs cattering  measurements.  TaSi2  is  equally  stable 
as  compared  to  the  Ta-silicides ,  however  there  have  been  no  measurements 
of  silicide  formation  in  evaporated  Ta  films  on  Si.  For  the  Ti/Si  and  Zr/Si 
systems  the  monosilicides ,  TiSi  and  ZrSi,  have  lower  free  energies  than 
the  disilicides  TiSi2  and  ZrSi,,.  This  suggests  that  before  the  formation 
of  the  disilicide  an.  intermediate  stage  might  be  observed.  There  is  no 
direct  evidence  for  this  although  our  preliminary  results33  on  Ti/Si  hint 
at  the  presence  of  an  intermediate  phase  before  Tisi2  is  formed,  in  the 
Hf/Si  system,  which  should  be  analogous  to  Zr/Si,  HfSi  does  form  and  is 
stable  up  to  temperatures  of  $  750«C  where  HfSi2  forms.  In  spite  of  these 
indications  we  believe  that  there  is  not  yet  sufficient  evidence  to  make  a 
strong  correlation  with  formation  energy.  This  situation  might  be  improved 
if  measurements  were  made  on  the  Ta/Si  system  and  further  measurements  were 
carried  out  in  the  Mo/Si  system. 
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From  a  conventional  standpoint,  the  stability  of  a  phase  is  indicated 

by  its  melting  point.  This  may  be  misleading  and  should  not  be  used  as 

a  guide  to  predict  which  forms  first.  It  is  interesting  to  note,  however, 

from  Table  I  that  a  value  of  half  the  melting  point  in  °K  is  rough  guide 

to  the  minimum  temperature  at  which  silicide  formation  has  been  observed. 

This  may  be  too  high  an  estimate  because  HfSi ,  Pd2Si ,  WSi2  and  PtSi  form 

at  about  one-third  the  melting  point.  In  fact  the  one-third  relation  may 

be  the  rule  rather  than  uhe  exception  because  oxide  layers  at  the  interface 

will  raise  the  formation  temperature.  For  example,  the  published  NiSi 

formation  temperature  of  600 °C  is  higher  than  that  of  300°C  recently 

34 

obtained  in  a  separate  investigation  .  We  suspect  that  the  presence  of 
the  native  oxide  layer  on  Si  might  be  responsible  for  the  high  temperature 
found  for  the  Mo/Si  system. 

We  can  consider  other  aspects  such  as  epitaxial  growth  that  might 

12 

effect  stability.  Hutchins  and  Shepala  showed  for  transformation  of 
Pd2Si  to  PdSi  that  higher  temperatures  are  required  on  <111>  oriented  Si 
than  on  <100>.  On  <111>  oriented  Si,  Pc^Si  shows  the  strongest  orientation. 
While  Pd„Si  is  very  stable  { 200oC-700oC)  ,  Pt  Si  is  not  and  transforms  to 
PtSi  quite  easily.  One  reason  might  be  that  Pd2Si  is  epitaxial  and  has 
a  higher  melting  point  than  PdSi  while  Pt^Si  does  not  show  good  epitaxy 
and  has  a  lower  melting  point  than  PtSi. 

For  Pt2Si  one  finds  the  low -temperature ,  (tetragonal)  phase  in  these 
silicides  rather  than  the  high  temperature  (hexagonal)  phase;  the  polymorphic 
transition  temperature  is  at  about  700°C.  We  attempted  to  form  this 
hexagonal  phase  (which  should  epitaxial  on  <111>  Si)  by  heating  a  sample 
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to  750 °C  for  30  min. 


Backs  cat tering  data^  showed  that  PtSi  (3%  Si  rich) 
rather  than  Pt2Si  was  formed.  Consequently  this  attempt  to  form  a  more 
stable  phase  of  Pt2Si  was  not  successful. 

C.  Kinetics 

From  the  standpoint  of  kinetics,  the  important  factors  are  the 
mechanisms  for  transporting  metal  or  Si  across  the  silicide  layer,  the 
activation  energy  and  the  identity  of  the  diffusing  species.  Most 
silicides  tend  to  react  at  temperatures  above  500-600°C  except  for  the 
Pd/Si  and  Pt/Si  systems  which  react  at  temperatures  a3  low  as  200 °C.  The 
activation  energy  measured  for  the  growth  process  for  these  latter  silicides 
is  between  1.1  and  1.5  eV.  This  is  remarkably  low  and  opens  the  question  of 
grain  boundary  growth.  We  believe  however  that  bulk  diffusion  dominates. 

The  structure  of  Pd2Si  is  very  °Pen  30(3  interdiffusion  is  expected  to  be 
fast.  Further  the  same  activation  energy  was  found  from  100°  to  700  C 
indicating  that  the  same  diffusion  mechanism  was  involved.  Also  the  same 
growth  rate  is  found  for  Pd2Si  on  different  oriented  Si  samples  where  the 
amount  of  epitaxy  is  different.  In  cases  where  there  is  less  epitaxy,  one 
anticipates  mors  large-angle  grain  boundaries  and  increased  grain  boundary 
diffusion.  All  these  findings  support  the  dominance  of  bulk  diffusion 
even  at  temperatures  of  200°C. 

The  other  silicides  react  -t  higher  temperatures  and  we  would  expect 
lattice  diffusion  to  be  dom.nant.  The  activation  energy  for  diffusion 
has  only  been  measured  in  a  few  cases  but  seems  to  be  above  2  eV  (2.5  eV 
for  HfSi ,  2.9  for  VSi2  and  2.7  for  WSi2) .  These  values  are  too  high  for 
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grain  boundary  diffusion  and  suggest  bulk  effects.  One  would  expect  to 
see  grain  boundary  diffusion  at  lower  temperatures.  This  could  be  tested 
by  prolonged  heating  of  a  sample  at  temperatures  below  the  reaction  tempera¬ 
tures  and  looking  for  Si  accumulation  at  the  free  surface.  This  interface 
accumulation  is  often  seen  in  backscattering  studies  of  metal-metal  inter¬ 
actions  at  low  temperatures.  However,  no  systematic  studies  have  been  made 
on  metal/Si  systems. 

The  growth  of  most  silicides  (HfSi,  Pd  Si,  Pt  Si,  PtSi,  TiSi  ,  WSi 

‘  2  2  2 

and  VSi2)  have  been  found  to  follow  an  (time) 1/2  relationship  characteristic 

of  a  diffusion  dominated  process.  In  two  cases,  CrSi  and  MoSi  a  linear 

2  2 
c. 

growth  has  been  found.  These  systems  have  not  been  studied  in  detail  and 
we  cannot  speculate  on  the  reasons  for  the  linear  growth  rate.  In  fact, 
the  growth  kinetics  of  all  the  high  temperature  silicides  has  not  been 
followed  in  as  much  detail  as  has  been  done  for  Pd/Si  and  Pt/Si.  There  is 

a  need  for  careful  kinetic  studies  on  samples  with  oxide  free  interfaces 
between  metal  and  Si. 

There  is  one  other  notable  gap  in  all  the  data.  With  the  exception 
of  HfSi ,  there  is  no  evidence  to  indicate  whether  Si  or  metal  is  the  diffusing 
species  through  the  silicide  layer.  This  data  could  be  obtained  by  implanting 
rare  gas  atoms  into  the  silicide  and  using  backscattering  techniques  to  determine 
whether  the  gas  atoms  moves  toward  or  away  from  the  surface.  Similar  studies  have 
been  made  with  anodic  oxidation.35  For  the  Hf/Si  system,  the  argon  that 
was  present  as  an  impurity  introduced  during  deposition  served  as  a  marker.7 
The  data  indicated  the  Si  was  the  diffusing  species.  There  are  practical 
implications.  If  si  is  the  diffusing  species,  it  will  leave  vacancies  near 
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voids,  it  leads 


the  Si/aiUcide  interface.  If  these  vacancies  condense  into 
to  easy  peeling  of  the  silicide  under  stress. 


D.  Comparison  with  M/Si02  Reactions 

It  has  recently  been  found  the  silicido 

slAlclde  formation  can  occur  with  metals 

11,16 


in  contact  with  siC>2. 


The  behavior  of  V/Si02  and  V/Si  has  been 


discussed  above  in  Sec.  IIB.  In  the  case  of  Ti,  Nb,  and  V  the  silicide 
was  the  more  mere!  rich  phase.16  At  present  our  work  is  in  too  early  a 
stage  to  give  any  i„ei,ht  into  this  behavior.  One  inflate  question  is 
to  determine  rf  the  Si  rich  phase  forms  when  the  Si02  layer  is  consumed. 


V  SUMMARY 

This  review  was  intended  as  a  survey  of  metal-silicide  formation 
studies  where  both  backsc.ttaring  and  x-ray  diffraction  measurements  have 
been  utilized.  These  techniques  complement  each  other  as  backscattering 
data  gives  composition  as  a  function  of  depth  and  hence  can  be  used  to 
determine  growth  kinetics  while  x-ray  diffraction  gives  identification 
Of  phases  and  structural  information. 

A  survey  of  silicide  formation  for  metal  films  deposited  on  Si  indi¬ 
cates  that  the  typical  silicide  that  is  formed  is  the  phase  that  is  most 
rich  in  Si.  ais  is  the  end  point  one  expects  from  thermodynamic  equilibrium 

considerations.  I"  seme  cases  intermediate  phases,  have  been  found.  On 

the  basis  of  present  data  we  have  not  been  able  to  find  general  guidelines 
that  can  be  used  to  predict  when  sud,  an  intermediate  phase  win  occur, 
one  striking  observation  is  that  metal  rich  silicides  are  formed  when  the 
metal  layer  is  deposited  on  Si02  rather  than  on  Si. 
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The  formation  temperature  on  the  silicides  is  generally  between  one- 
third  to  one-half  the  silicide  melting  point  in  °K.  This  shews  that  solid- 
solid  rather  than  solid-l.iquidus  reactions  are  involved.  The  growth 
kinetics  of  the  silicide  with  two  exceptions  have  been  found  to  follow  a 
(time) 1/2  dependence  typical  of  diffusion  limited  reactions.  We  believe 
that  bulk  rather  than  grain  boundary  diffusion  is  involved.  There  have 
not  been  studies  to  determine  whether  the  metal  or  the  Si  is  the  diffusing 
species  except  for  the  Hf/Si  case  where  Si  appeared  to  migrate. 

Our  survey  indicates  that  only  the  most  general  features  of  silicide 
formation  have  been  studied.  There  are  a  number  of  areas  that  require 
further  investigation.  Until  these  are  carried  out  it  is  difficult  to 
present  an  overall  picture  of  the  factors  which  control  silicide  formation 
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SILICIDE  FORMATION 


Backscatterlnn _  obs 


t'etal+ 

Phases 

Structure 

Phase 

T  ,  *C 

O.'S 

Kinetics 

T  -C 

melt 

T  *K 

melt 

Ref 

Ti5Si3 

hexagonal 

— 

2120 

T1  (3) 

TiSl 

orthohombic 

— 

1700 

TISlj 

(C49) 

TlSl2 

600 

tU2 

1540 

% 

-v-0.5 

6 

Zr-Sl 

ZrSi 

tetragonal 

-- 

.  2110 

Zr  (7) 

ortho 

— 

2095 

15 

ZfSij 

"  (C49) 

ZrSij 

700 

-  * 

1520 

So.5 

Hf,Si 

tetra 

2430 

Hf  (5) 

2 

HfSi 

ortho  or  hexa 

Hf  Si 

550 

tl/2 

2200 

% 

Ml. 3 

(FeB) 

7,8 

HfSlj 

ortho  (C49) 

HfSij 

750 

1900 

- 

M).5 

V  SI 

cubic  (9-W) 

2070 

V  (3) 

V.Si, 

tetra  (D8m) 

2150 

5  3 
VSi2 

hexa  (C40) 

vsi2 

600 

tl/2 

1750 

M).5 

11,16 

Nb.Sl  hexa 

4 

Nb  (4)  NbjSl  cubic  (Cu.jAu) 

NbSlj  hexa  (C40)  NbSij  650  ^0.5  30 


Ta.Si 

tetra  (Al,Cu) 

2460 

Ts  (3) 

2200 

TaSlj 

hexa  (C40) 

Cr3Si 

cubic  (fl-W) 

1730 

Cr  (5) 

CrSl 

cubic  (FeSi) 

1600 

6,15 

CrSij 

hexa  (C40) 

CrSi2 

650 

t 

1550 

Ml. 4 

Mo.Si 

cubic  (0-W) 

2120 

Mo  (3) 

tetra  (cilb) 

MoSij 

MoS12 

1200 

t 

2050 

■vO.6 

6 

W  St, 

tetra  . 

2350 

W  (2) 

J  L 

WS12 

tetra  (Cll) 

wsi2 

650 

t1/2 

2165 

So.35 

17 

NljSl 

cubic  (CUjAu) 

1165 

Ni  (6) 

NISI 

ortho  (MnP) 

NISI 

600 

992 

_ 

M).7 

NiSlj 

cubic  (CaFj 

NiSlj 

800 

993 

^0. 85 

15 

Pd  SI 

ortho  (Fe.C)  ' 

t1/2 

960 

Pd  (3) 

Pd,Sl 

hexa  (Fe^F) 

Pd, SI 

200 

1330 

So.  35 

3,6,9,10 

12,13,14 

PdSl 

ortho  (MnP) 

PdSl 

735 

1090 

So. 75 

14 

Pt  SI 

mono 

t1/2 

870 

Pt  (5) 

pt2si 

PtSl 

hexa  (Fe?P) 
tetra  (AT.Cu) 
ortho  (MnP) 

Pt2Si 

PtSl 

200 

200 

1100 

1229 

M>.4 

% 

M>.35 

4,5 

+Tha  nuaber  In  parentheses  Indicated  the  number  of  phases  fotnvl  In  bulk  samples.  (18-20) 
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Schematic  diagram  of  the  sample  chamber  used  in  backscattering 
measurements . 
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Backscattering  spectra  for  2  MeV  He  ions  incident  on  Si  sample 
covered  withc3.l00A  of  V  and  heat  treated  at  600 °C.  Data  from 


H.  Kravitle 
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Figure  3.  B.ckscatterlng  spectra  for  2  MeV  4He  ion.  incident  on  a  ”^i- 
layer  antireflecting  coating  on  quartz  substrate.  There  are 
layers  of  ThF  separated  by  layers  of  ZnS  (sample  supplied  by 

R.  Honig,  RCA). 


3200  &  V/Si  800°C  15  MIN 


Figure  5 


X-ray  diffraction  pattern  of  V/Si  heat  treated  at  800°C  for  13 
minutes  in  an  oxygen  free  furnace.  The  vanadium  film  had  a 
thickness  of  3200A.  Reflection  of  VSi_  are  indexed. 
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Figure  6.  Backs cattering  spectra  for  2.3  MeV  4He  ions  incident  on  a  Si 

sample  with  a  thermally  grown  layer  of  Sf°2.  °°^oc  fot  60^dn  in 
evaporated  layer  of  V.  A)  heat  treatment  at  400  C  for  60  Bin.  in 

0  ambient,  and  B)  treatment  at  800 °C  for  30  min.  in  v 

fyrtm  M  Krm'itlfi 


X-ray  diffraction  pattern  of  V/SiO  heat  treated  at  400»c  for  4 
hours  in  an  atmosphere  of  0  +  H  O?  The  vanadium  film  had  a 
thickness  of  3000A.  Reflection  of  v„0_  are  indexed. 
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SUICIDE  FORMATION 


ZrSI,  (700) 


NbSI-  (650) 


MoSI,  (1200?) 


HfSI  (550) 

Pt2S1  (200° ) 

IfSI,  (750) 

WS12  (650°) 

PtSI  (300°) 

Figure  8.  A  conpariaon  of  silicidc  formation  n  thin  metal  films  deposited 
on  Si.  the  numbers  in  parentheses  indicate  the  temperature  at 
which  silicide  formation  was  observed  by  backs  cat tering  techniques 
The  asterixes  indicate  silicides  that  been  identified  by 

diffraction  techniques. 
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ABSTRACT 

The  changes  in  vacuum-evaporated  films  of  Ti,  V  and  Nb  on  Si  and 
Si02  substrates  after  thermal  anneals  are  investigated  by  backscatterinq 
and  by  x-ray  spectrometr  Backscatter ing  analysis  provides  the  relative 

atomic  composition  as  a  function  of  depth  with  high  sensitivity.  Glancing 
x-ray  spectrometry  detects  the  chemical  composition  with  high  specificity. 
Combined,  the  two  methods  create  a  specific  picture  of  the  transformations 
induced  in  the  films  by  the  thermal  treatment.  Generally  the  reaction  on 
a  pure  Si  substrate  produces  a  si-rich  silicide  and  on  a  Si02  a  metal  rich 
silicide  layer  as  an  intermediate  layer  largely  free  of  oxygen.  The 
oxygen  originally  bound  to  the  Si  in  the  Si02  is  transferred  to  the 
remaining  metal  layer.  Residual  oxygen  in  the  metal  film  and  metal  oxides 
on  the  metal  film  influence  the  silicide  formation. 
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INTRODUCTION 


Transition  metals,  and  Ti  in  particular,  are  used  extensively  in 
metallization  schemes  for  integrated  circuits  and  solar  cells.  Good 
adhesion  and  uniformity  of  the  resulting  contact  are  main  reasons  for 
their  frequent  application.  The  latter  aspect  has  recently  been  inves¬ 
tigated  in  connection  with  the  Al-Ti  metallization  scheme. (1)  The 
possibility  to  form  superconducting  silicides  of  transition  metals  by 
thin  film  reactions  with  the  substrate  has  prompted  Tu ,  et  al .  to  inves¬ 
tigate  the  behavior  of  V  on  Si  and  Si02  substrates .  ( 2 1  The  strong 
affinity  of  Ti  and  V  for  oxygen  is  believed  to  be  the  reason  why  these 
metals  make  good  blocking  contacts  to  p-type  semiconductor  oxides  and 
good  ohmic  contacts  to  n-type  oxides.'31  Niobium  (Columbium)  maintains 
high  mechanical  strengths  at  elevated  temperatures,  hut  oxidizes  reaaily 
unless  protected  by  a  coating.  The  disilicide  constitutes  such  a  pro¬ 
tective  layer,  based  on  the  ability  of  the  coating  to  generate  a 
silica-  ed  glassy  oxide  as  protection  against  the  atmosphere. (4) 

Preliminary  data  indicate  that  the  heat  treatment  of  these  metal 
films  on  Si  results  in  the  formation  of  silicides,  but  that  the  same 
metal  films  deposited  on  Si02  forms  both  silicides  and  oxides,  generally 
separated  in  distinct  layers.'51  we  present  here  a  systematic  study 
of  samples  annealed  in  both  vacuum  and  in  dry  ox-.en.  Comparative 
investigations  such  as  these  should  lead  to  insight  in  the  dominant 
processes  at  work  in  these  thin  films,  and  to  a  better  understanding 
of  their  applications. 
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II.  METHODS  AND  ANALYSES 


A.  Backscattering  Spectrometry 

In  this  study,  the  depth  and  mass  perception  of  backscattering 
analysis  are  used  to  determine  the  relative  atomic  composition  of  the 
films.  The  backscattering  technique  and  the  brief,  a  monoene rgetic  and 
well  collimated  Learn  of  He  ions  with  energy  ranging  from  2.0  to  2.4 
MeV  impinges  r erpendicularly  onto  the  sample.  A  small  fraction  of 
these  ions  is  scattered  back  into  a  surface  barrier  Si  detector  which 
is  mounted  at  an  angle  of  168  degrees  against  the  incident  beam.  The 
detector  signals  are  amplified,  shaped,  and  recorded  in  a  multi-channel 
analyzer.  The  resulting  energy  spectrum  (counts  per  energy  interval 

versus  energy)  furnishes  the  information  of  atomic  composition  versus 
depth . 

As  an  example.  Fig.  la  shows  an  energy  spectrum  of  2  MeV  4He+  back- 

scattered  from  2100A  of  vandium  deposited  on  2700A  of  thermally  grown 

Si02  on  a  Si  substrate.  Helium  ions  backscattered  from  oxygen,  silicon 

and  vanadium  atoms  in  the  sample  have  energies  which  depend  on  the 

mass  of  the  scattering  atom  and  its  location  in  depth.  After  the  same 

sample  is  annealed  in  vacuum  at  800*C  for  two  hours,  the  backscattering 

analysis  is  repeated  and  the  energy  spectrum  given  in  Fiq.  lb  is 

obtained.  For  ease  of  interpretation,  various  areas  under  the  energy 

spectra  which  correspond  to  different  layers  have  been  shaded  (see 

also  insert  in  :he  figures).  The  two  spectra  clearly  differ,  as  can 

be  seen  in  Fig.  lb  where  for  comparison  the  spectrum  of  Fig.  la  is 

replotted  as  a  dashed  line  The  formation  of  two  distinct  layers 

(V  Si  and  VO  )  and  the  remaining  part  of  the  SiO 
*  2 


is  readily  distinguished. 


ENERGY (MeV) 

Figure  1.  Spectrum  of  2  MeV  lle+  backscattered  from  a  V  film  on  an  SiC>2 
layer  on  Si  as  evaporated  (top),  and  after  heat  treatment 
in  vacuum  at  800°C  for  2  hours  (bottom) . 

One  can  convert  a  backscattering  energy  spectrum  into  an  atomic 
concentration  profile  with  the  energy  loss  factor  ( 5  (  ■  AE/Ax  which 
relates  a  change  in  AE  in  backscatter  i'.y  energy  to  a  change  Ax  in  the 
depth  of  the  sample.  ^ ^  The  (S)-factor  depends,  among  others,  on  the 
energy  of  th  incident  particle  and  the  composition  of  the  layer.  The 
energy  dependence  is  slow  and  for  thin  films,  ( S ]  can  be  assumed  to  be 
a  constant.  To  a  good  approximation,  the  ratios  of  the  [Sj  factors  for 
the  various  elements  in  a  film  is  also  independent  of  composition.  By 
measuring  the  ratio  of  the  signal  height  generated  in  a  spectrum  by  the 


elements  of  a  fil..  one  therefore  obtains  the  elements!  concen- 

tios  Of  this  unknown  fil-.  *  correction  factor  must  be 
tration  ratios  Rutherford  scattering 

„  for  the  different  (but  known) 
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can  readily  identify  chemical  comf«sltlon.  Bec.09e 

.  .fll,e  to  backscattering  analysis, 
silent  complementary  techn  q 
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incident  angle.  Per  example,  an  incident  angle 
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th.  P,th  length  of  th.  .....  in  the  to  .bout  12  time.  it.  thick- 

.....  th.  focusing  Seemann-Bohl in  arrangement  th.  sp.ci..n  is  pl.o.a 

„  th.  circumf.t ence  of  th.  diffraction  citcl..  Th.  .ngl.  of  incid.no. 

is  fixed  ..a  can  b.  ..a.  ..  —H  ■■  *  °"  ^ 

principle ,  ».  r.d.r  ..a  ..  ..  Bert,  have  d.slgn.d  .nd  built  .  ^inches 

diameter  x-ray  diff r.cto.e t.r  .pecifi-HV  ^  thin  film  studies  . 

Th.  dif f ractom. t.r  ..ploy.  .  plyrolytlc  gr.phit.  monochrom.t.r  crystal 

to  obtain  a  high  intensity  monochromatic  cuka  radiation,  which  is 

incident  upon  th.  specimen,  at  an  angle  of  6.4*. 

This  diffractometer  possesses  sufficient  sensitivity  to  resolve 
the  first  five  diffraction  peaks  of  a  copper  film  of  only  150h  thic„..ess 
using  a  scan  in  et.PS  of  0.2-  ,40,  and  a  counting  time  of  75  sec/step. 

The  broadening  exhibited  by  the  diffraction  profile  is  consistent  with 
a  particle  six.  of  less  than  lool.  The  diffractometer  hold,  great 
promise  in  the  investigation  of  the  kinetic,  of  reaction  in  thin  films, 
hecaus.  th.  reaction  product,  can  be  identified  b,  their  x-ray  reflect¬ 
ion.  very  early  in  the  diffusion  cycle,  and  the  rat.  of  change  can  be 
measured  by  changes  of  their  peak  intensities."”  The  capability  of 
the  diffractometer  to  determine  th.  lattice  parameter  very  precisely 
.1st  allows  the  measurement  of  stress  in  a  thin  film.  The  measured 
lattice  parameter  of  a  well  annealed  hi  Powder  ...Pie  -«  1* 
agreement  with  the  value,  a  -  ,26-C,  given  in  the  hSTh  file, 

using  the  powder  sample  as  the  reference,  the  stress  state  of  a  number 
loooi  hi  films  on  glass  substrates  was  determined  and  found  to  ^ 
X.lt„e  with  th.  results  obtained  by  the  bending  contilev.r  technigue 
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V  on  Si02  (90CTC  I  HR.) 
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Figure  2.  X-ray  diffraction  pattern  of  the  sample  of  Fig.  1  after  heat 
treatment.  Reflections  of  V^Si,  V,.Si3  end  v2°s  are  observed* 

We  have  used  the  diffractometer  to  study  the  reaction  products  of 

our  transition  metal  films.  The  x-ray  analyses  were  performed  partly 

on  samples  identical  with  those  measured  by  backscattering,  partly  on 

samples  identical  with  those  measured  by  backscatte r  ir.g ,  partly  on 

samples  prepared  and  annealed  under  closely  similar  conditions.  For 

2 

best  results,  samples  of  several  cm  area  are  required  for  x-ray  measure 
ments ,  which  is  typically  10  times  the  size  needed  for  backscatterinq 

analysis.  Figure  2  shows  a  typical  x-ray  spectrum  for  a  sample  of 

•  • 

2000A  of  V  on  an  oxidized  Si  wafer  (5000A  of  SiO^)  after  anneal  at 

800  *C  for  30  minutes  in  vacuum.  The  spectrum  was  obtained  by  scanning 
the  sample  at  steps  of  0.15*  (40)  increment  and  with  a  counting  time 
of  30  sec/step.  The  peaks  have  been  indexed  as  reflection  of  VSi^, 


V c S i  and  V_0  . 


C.  Sample  Preparation 

For  the  experiments  on  Si  substrates,  commercially  polished  wafers 
of  either  <100>  or  <111>  orientation  and  of  usual  n-  and  p-type  doping 
levels  were  employed.  Orientation  and  dopinq  have  no  detectable 
influence  on  the  results.  Immediately  before  vacuum  evaporation  of 
the  metal  film,  the  wafers  were  etched  in  a  solution  of  2  HNO^l  HF:2 
acetic  acid  and  rinsed  in  deionized  water.  The  SiO.,  substrates  were 
obtained  by  thermally  oxidizing  polished  Si  wafers  at  1100‘C  in  a 
atmosphere  of  wet  oxygen.  All  metal  films  were  deposited  by  vacuum 

e vaporat ion . 

III.  REACTIONS  WITH  Si  AND  SiOj  SUBSTRATES 
A.  Titanium 

The  result  of  the  reaction  of  Ti  with  a  Si  substrate  can  be  seen 
in  Fig.  3.  The  zero  in  the  depth  scale  is  placed  at  an  arbitrary 
reference  point  within  the  substrate,  and  positions  between  this  point 
and  the  surface  are  indicated  with  thicknesses  greater  than  zero. 

Figure  3  (top)  shows  2700A  of  Ti  as  evaporated  on  a  Si  substrate.  The 
bottom  part  of  this  figure  shows  this  sample  after  heat  tJ eatment  at 
600*C  for  20  min  in  a  vacuum  of  bettei  than  10  Torr.  The  uppermost 
layer  now  contains  oxygen  which  wall  absorbed  by  the  Ti  layer  during  the 
storage  in  air  and  also  during  annftaling.  The  spread  of  the  points 
exhibits  the  poor  sensitivity  of  the  backscattering  method  for  elements 
of  lower  mass  than  that  of  the  substrate.  Below  the  Ti  on  the  surface, 
an  intermediate  layer  has  been  formed  during  anneal.  The  ratio  of  Ti 


to  Si  in  that  layer  can  k 
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The  chemical  structure  has  be  ,  18  t0  **  1:2 
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The  interaction  of  a  vacuum-evaporated  Ti  film  with  an  amorphous 
Si02  substrate  is  more  complicated,  since  the  metal  reacts  with  a 
compound.  A  good  example  of  the  complexity  of  this  system  is  shown  in 
Fig.  4.  The  top  figure  gives  the  concentration  profile  of  a  sample 
after  deposition  of  the  metal  film,  and  the  figure  below  shows  what 
happens  to  that  sample  after  heat  treatment  at  800‘C  for  2  hours  in 
vacuum.  The  surface  layer,  which  was  pure  Ti  before  annealing,  contains 
about  50%  oxygen  after  annealing.  The  layer  beneath  consist  mostly  of 
Si  and  Ti  in  a  ratio  Si:Ti  -  0.62.  A  small  amount  of  oxygen  is  present 
also.  Its  concentration  is  difficult  to  evaluate  by  the  backscatter ing 
method.  Within  an  accuracy  of  about  10%  the  spectrum  tells  that  the 
amount  of  oxygen  in  the  surface  layer  is  equal  to  that  originally 
present  in  the  Si02  which  reacted  with  Ti  to  form  a  silicide.  X-ray 
analysis  could  not  determine  the  composition  of  this  Ti-0  layer. 


B.  Vanadium 

Far  clearer  is  the  interaction  of  V  films  with  Si  and  Si02.  In 
Fig.  5  the  reaction  of  V  with  pure  Si  can  be  seen.  The  intermediate 
layer  which  grows  between  the  Si  substrate  and  the  V  film  shows  two 
well  delineated  interfaces  in  the  backscattering  spectrum.  There  is 
no  detectable  oxidation  of  V  from  the  air  at  room  temperature  or  during 
annealing  in  our  vacuum  system.  The  backscattering  spectrum  indicates 
a  ratio  of  V:Si  =  1.2  in  the  intermediate  layer.  The  x-ray  analysis 
substantiates  the  result.  Only  one  compound  VSi2  is  formed  over  the 
whole  temperature  range  from  *00-C  to  1000-C  In  addition,  the 
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relative  intensity  ratios  of  tie  x-ray  diffraction  pattern  indicate 
no  preferred  orientation  in  the  silicide. 

At  emperatures  of  700#C  and  above  V  films  react  with  Si02,  as  can 
be  seen  in  Fig.  6.  The  various  layers  observed  after  the  reaction 
are  .eery  well  resolved  in  the  backscatter  ing  spectra  (see  Fig.  l). 
According  to  this  spectrum  ".he  surface  layer  after  reaction  consists 
of  V  and  0  in  the  ratio  of  V:0  about  1:1.  According  to  the  x-ray  data, 
however,  tho  layer  contains  V.,0,.  and  V,.Si3.  No  measurable  amount  of 
Si  has  been  found  in  this  layer  by  Rutherford  backscattering  analysis. 
The  intermediate  layer  contains  only  V  and  Si.  The  composition  is 
V.Si  ■  3:1  from  back r cat ter ing  data  and  is  identified  as  pure  V  Si  by 
x-ray  analysis.  The  conservation  of  the  total  oxygen  content  in  the 
system  is  verified  within  an  accuracy  better  than  10%  by  comparing 
the  total  oxygen  amount  in  layers  before  and  after  heat  treatment,  or 
by  measuring  the  position  of  the  rear  edge  of  the  oxyqen  signal  of  the 
Si02  layer.  This  edge  must  move  if  any  substance  evaporates  on  or  off 
the  surface,  or  diffuses  in  or  out  of  the  system.  This  indicates  also 
that  no  measurable  amount  of  the  substrate  (Si)  diffuses  through  the 
Si02  layer  into  the  substrate. 

C.  Niobium  (Columbium) 

The  reaction  temperature  for  Nb  films  on  Si  and  on  SiO  prepared 
at  Caltech  are  higher  than  700*  and  900*C  respectively.  Unfortunately, 
several  experimental  flaws  affect  the  results.  The  evaporated  layers 
always  contain  some  oxygen.  During  heat  treatment  the  oxygen  content 
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Figure  5.  (left)  Concentration  profile  of  a  V  film  on  Si  as  evaporated 
(top),  and  after  heat  treatment  in  vacuum  at  600“C  for  20 
minutes  (bottom) . 

Figure  6.  (right)  Concentration  profile  of  a  V  film  as  evaporated  on 
an  oxidized  Si  substrate  (top),  and  after  heat  treatment  in 
vacuum  at  800°C  for  30  minutes  (bottom) . 

increases  further.  From  the  best  evaporations,  we  determine  for  both 

substrates  a  ratio  for  Nb:Si  of  about  5:3  for  the  intermediate  layers. 

Since  we  did  not  succeed  in  preparing  or  maintaining  oxygen-free  Nb 

layers,  the  top  layer  of  Nb  has  always  a  noticeable  oxygen  content. 


IV.  OXIDATION  OF  METAL  FILMS  ON  Si02 
The  oxidation  products  of  the  metal  layers  in  an  oxygen  atmosphere 


have  also  been  studied. 


Tha  results  can  be  seen  in  Figs.  7-9. 


Ti  is  known  as  a  very  reactive  material  able  to  incorporate  a 
remarkable  amount  of  oxygen  even  at  room  temperature.  We  have  observed 
that  heat  treatment  in  vacuum  tends  to  homogeneously  distribute  oxygen 
throughout  the  Ti  film,  indicating  that  oxygen  is  very  mobile.  Heat 
treatment  in  dry  oxygen  atmosphere  at  400*C  for  2  hours  show  a  total 
oxygen  concentration  of  about  50%  which,  however,  is  not  uniform  dis¬ 
tributed  in  the  layer.  Hut  at  600*C,  a  rapid  increase  of  the  oxygen 
concentration  can  be  seen  (Fig.  7,  lower  part).  The  layer  can  be 
divided  into  two  regions;  the  top  layer  rhows  an  oxygen  to  Ti  ratio 
of  about  3:2  and  the  adjoining  region  below,  which  does  not  reach 
visibly  with  the  SiC>2  layer  at  that  temperature,  has  a  ratio  0:Ti  “  1:1. 

Vanadium  seems  to  form  a  compound  at  relatively  low  oxidation 
temperatures.  In  Fig.  0,  the  composition  of  the  sample  can  be  seen 
before  (top)  and  after  (bottom)  heat  treatment  at  400°C  for  30  minutes 
in  dry  oxygen.  The  transition  between  the  oxide  layer  and  the  V  metal 
layer  is  sharp.  The  ratio  of  the  oxidized  layer  is  V:0  -  2.5,  and  is 

identified  by  x-ray  analysis  as  v2°5* 

This  compound,  which  is  the  same  as  that  found  by  x-ray  diffraction 
in  the  product  of  the  SiO^-V  reaction,  has  as  low  melting  point 
(660*0  To  evaluate  the  effects  associated  with  this  melting  point, 

the  following  experiment  has  been  performed.  A  V2C>5  layer,  formed  by 
complete  oxidation  of  a  V  film  on  Si02  at  600*C  in  dry  oxygen,  was 
further  heat  treated  at  700«C  in  oxygen.  The  layer  melts  and  forms 
drops  on  the  surface.  The  same  V2C>5  layer  annexed  below  600*C  in 
vacuum  shows  no  visible  change.  But  during  an  anneal  at  640*C  in 
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Figure  7.  (left)  Concentration  profile  of  a  Ti  film  as  evaporated  on 
an  oxidized  Si  substrate  (top),  and  after  heat  treatment  in 
dry  oxygen  at  600*C  for  one  hour  (bottom) . 

Figure  8.  (right)  Concentration  profile  of  a  V  film  as  evaporated  on 
an  oxidized  Si  substrata  (top),  and  after  heat  treatment  in 
dry  oxygen  at  400»C  for  30  minutes  (bottom). 

vacuum  for  several  hours,  most  of  the  VjOg  layer  vaporizes.  The  behavior 

of  an  only  partly  oxidized  V  film  on  Si02  as  shown  in  Fig.  8  is  different 

Up  to  600 *C  in  vacuu  1  no  change  can  be  seen.  But  at  temperatures  higher 

than  700*c  the  top  layer  changes  composition .  Tho  V:0  ratio  after 

annealing  for  several  hours  is  abo«t  1:1  and  below  this  thick  V-0  layer 

is  a  thin  V3Si  layer  from  the  reaction  of  Si02  with  the  excess  V.  The 

V  and  oxygen  content  in  conserved  within  the  accuracy  of  the  measurement. 

The  same  sample  ha*  been  annealed  again  in  a  dry  oxygen  atmosphere  at 


400 *C  and  the  top  layer  oxidizes  until  a  ratio  of  V:0  •  2:5  is  reached. 


Figure  9.  Concentration  profile  of  a  Nh  *41. 

oxidized  Si  substrate  < t'oj,  ?  ^ e’r TtlV Vrll t^n l\  d 

oxygen  at  400*C  for  one  our  (bottom).  treatment  in  dry 

annealing  of  a  Nb  film  on  SiC>2  in  oxygen  at  400*C,  a  layer 
of  the  composition  Nb:0  -  2:5  can  found  (Fig.  9). 


v.  discussion  and  conclusion 


Th.  bin.ry  P,...  dl„r...  . . .  ,ho-  th>t  a 

V*tUty  #*  ,1Ue,4“  “»  «*«•  *-  .v.t.m  four  co.poll„d, 

Ti5Sij,  Tis  1 ,  Tisip,  for  tb.  sl-v  system  three  compound, 

,v3Si.  v5sis,  vsi2)  ,„d  for  the  Si-Nb  system  four  compound,  (»b  Si, 

»b  Si,  Nb  Si  ,  NbSi  )  ,re  reported ,  110  • 11 1  * 

3  J  2  F  No  Prediction  can  be  made 

whet  compound  .in  be  formed  firet,  if  .„y,  by  i„t.r.cti„„  of  the 
elements  or  by  int.r.ction  of  different  co.poundo  et  t.mp.r.ture.  b.lo. 
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A  summary  of  the  compounds  observed  in  the  reaction  products 
between  the  metal  and  the  Si  or  Si02  substrates,  and  with  oxygen  is 
given  in  the  Table  1.  We  list  the  ratio  of  the  elements  contained  in 
the  films  as  derived  from  backscatter ing  analysis.  Beneath  the  back- 
scattering  data,  the  chemical  compounds  are  listed  which  are  identified 
by  glancing  angle  x-ray  measurements.  An  approximate  temperature  is 
also  given  above  which  a  noticeable  reaction  rate  is  observed.  Cor¬ 
responding  data  of  x-ray  and  of  backscatter ing  analyses  were  all 
obtained  on  the  identical  sample,  except  where  marked  by  *.  The 
results  of  t».s  presert  expjriments  lead  to  the  following  conclusions: 

1.  Reaction  betweena  SiC>2  substrate  and  the  metal  films  requires 
temperatures  which  are  about  200«C  higher  than  those  needed 
for  the  reaction  of  a  Si  substrate  with  the  same  metal  film. 

2.  Sillcides  formed  between  metals  and  Si  are  generally  silicon 
rich  (e.g.  VS i 2  )  . 

3.  Silicides  formed  between  metals  and  Si02  are  generally  metal- 
rich  (e.g.  V3Si) . 

A.  When  metals  and  SiC>2  react  in  vacuum,  metal  silicides  and  metal 
oxides  are  both  formed,  but  in  separate  layers. 

5.  Metal  oxides  produced  by  reactions  with  SiC>2  have  less  oxygen 
than  the  metal  oxidized  formed  in  an  oxidizing  ambient. 

6.  The  total  amount  of  Si,  oxygen,  and  metal  present  in  a  sample 
does  not  change  when  the  reaction  takes  place  in  an  inert 
surrounding . 

7.  No  indication  of  reversibility  of  the  reaction  has  been 
observed . 

Questions  such  as  how  fast  the  reaction  proceeds,  what  terminates 
it,  and  what  happen  when  a  given  layer  is  completely  consumed  are 
difficult  to  answer .  Preliminary  experiments  indicate  that  some  of 
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the  complications  are  due  to  problems  of  the  cleanness  at  the  interface, 
and  oxygen  contamination  througout  the  layers. 
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